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ABSTRACT 

This  thesis  takes  five  FORTRAN  IV  programs  from  "Com- 
puter Programs  for  Computational  Assistance  in  the  Study 
of  Linear  Control  Theory"  by  Melsa  and  Jones  and  translates 
them  into  a  microcomputer  BASIC  language  to  run  on  an  inex- 
pensive microcomputer  system.   Three  of  the  five  programs  are 
state  variable  programs.   They  are  BASMAT  for  basic  matrix 
manipulation,  RTRESP  for  rational  time  response,  and  GTRESP 
for  graphical  time  response.   Two  are  transfer  function  pro- 
grams, FRESP  for  frequency  response  and  RTLOC  for  root  locus. 
A  user's  guide  and  example  are  included  for  each. 

A  final  example  is  used  to  demonstrate  the  utility  of 
the  two  transfer  function  programs  as  an  aid  to  direct 
digital  design  in  the  w' -plane. 
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I.   INTRODUCTION 

The  purpose  of  this  thesis  is  two  fold.   First,  it 
is  to  show  that  standard  computer  programs  useful  in  the  study 
of  linear  control  theory  may  be  adapted  to  run  on  an  inex- 
pensive home /microcomputer .   To  demonstrate  this,  five 
programs  were  chosen  from  [Ref.  1],  three  state  variable 
programs  and  two  transfer  function  programs.   The  three  state 
variable  programs  are  BASMAT ,  a  basic  matrix  manipulation 
program,  RTRESP,  a  rational  time  response  program,  and  GTRESP , 
a  graphical  time  response  program.   The  two  transfer  function 
programs  are  FRESP ,  a  frequency  response  program,  and  RTLOC , 
a  root  locus  program.   These  programs  as  they  appear  in  [Ref. 
1]  are  written  in  basic  FORTRAN  IV  language  to  be  run  on 
a  main  frame  computer  system  utilizing  standard  graphics 
subroutines.   These  programs  were  modified  and  rewritten  in 
a  microcomputer  BASIC  language  which  is  an  interpreted 
language.   Generally  these  programs  are  limited  to  systems 
of  10th  order.   It  is  felt  that  this  limitation  is  more  than 
acceptable  for  the  purposes  of  this  thesis.   Also  no  major 
effort  has  been  made  to  analyze  or  improve  the  efficiency  of 
the  numerical  methods  used.   An  effort  of  this  type  is 
advised  if  these  programs  are  to  be  modified  for  higher 
order  systems. 


The  second  part  of  this  thesis  investigates  the  methods 
and  relationships  involved  in  direct  digital  design  in  the 
w' -plane.   The  two  transfer  function  programs  adapted  to  run 
on  a  microcomputer  from  the  first  part  of  this  thesis,  FRESP 
and  RTLOC ,  are  used  to  aid  in  this  investigation. 

Section  II  summarizes  the  common  problems  and  consi- 
derations involved  in  the  translation  of  the  five  programs 
to  be  run  on  a  microcomputer  system.   Section  III  describes 
the  three  state  variable  programs,  BASMAT ,  RTRESP,  and 
GTRESP,  and  gives  an  example  of  their  use  and  output. 
Section  IV  is  similar  to  section  III  and  describes  the  two 
transfer  function  programs,  FRESP  and  RTLOC. 

Section  V  deals  with  the  w1 -plane.   Subsection  A  gives 
some  background  on  the  w' -plane  and  subsection  B  provides  a 
simple  example  using  the  two  transfer  function  programs  to 
compare  the  s  and  w1  planes.   Subsection  C  develops  templates 
of  some  constant  parameters  in  the  w' -plane.   Subsection  D 
ends  the  section  with  a  more  involved  example. 

Section  VI  provides  some  conclusions  and  recommendations. 


II.   TRANSLATION  CONSIDERATIONS 

The  five  programs  in  this  thesis,  BASMAT,  RTRESP,  GTRESP, 
FRESP,  and  RTLOC ,  were  translated  directly  from  programs  of 
the  same  name  in  [Ref.  1].   These  programs  were  written  in 
the  basic  FORTRAN  IV  computer  language,  a  compiled  language 
designed  to  run  on  a  main  frame  or  mini  computer  system.   In 
the  first  part  of  this  thesis  these  programs  are  translated 
into  a  microcomputer  BASIC  language  adapted  to  run  on  the 
microcomputer  system  described  in  Appendix  A. 

In  general  the  translation  from  FORTRAN  IV  to  the  micro- 
computer BASIC  used  posed  no  serious  problems.   Even  though 
the  BASIC  language  used  is  by  necessity  an  abridged  version 
of  the  BASIC  language  found  on  larger  more  expensive  systems, 
it  was  extensive  enough  to  provide  the  necessary  commands 
for  these  programs. 

The  input  portion  of  the  programs  were  rewritten  to  be 
interactive  for  convenient  keyboard  entry  of  problem  para- 
meters eliminating  input  formatting  errors.   The  programs 
requiring  more  extensive  input  were  given  the  added  capability 
of  saving  and  retrieving  a  problem  description  from  a  disk 
file  for  the  user's  convenience. 

Unlike  FORTRAN  IV  the  BASIC  used  does  not  have  complex 
math  capabilities.   Therefore,  mathematical  operations  on 
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complex  quantities  were  programed  separately  for  the  real  and 
imaginary  parts. 

Another  important  difference  in  the  languages  affecting 
translation  is  that  BASIC  has  no  provision  for  local  variables. 
All  variables  in  the  main  program  and  all  subroutines  are 
global.   This  caused  some  bookkeeping  problems  in  translating 
subroutines  to  prevent  undesired  side  effects.   Another  compli- 
cation was  that  FORTRAN  IV  considers  the  first  four  characters 
of  a  variable  name  for  identification  of  that  unique  variable. 
The  BASIC  used  considers  only  the  first  two.   This  further 
complicated  the  bookkeeping  of  variables  and  resulted  in 
variable  names  being  assigned  just  because  they  were  different 
from  the  rest  and  with  no  relation  to  the  quantity  represented. 

The  output  routines  were  written  to  conform  as  close  as 
possible  to  the  FORTRAN  IV  version.   Since  the  microcomputer 
system  used  is  limited  to  eighty  columns  on  eight  and  a  half 
inch  wide  paper,  provisions  were  made  within  the  program  to 
automatically  switch  to  a  condensed  character  font  when 
necessary  to  output  greater  than  eighty  characters  per  line. 

The  programs  requiring  graphical  output  created  some 
unique  problems.   First,  no  standard  library  subroutines  for 
graphics  as  used  in  the  FORTRAN  IV  programs  were  available 
on  the  microcomputer  system  used.   As  a  result  all  graphics 
routines  had  to  be  developed  for  the  system  used.   Also  the 
appropriate  variables  were  redimensioned  giving  consideration 
to  the  resolution  of  the  graphics  available  on  the 

10 


microcomputer  system  to  optimize  the  program  somewhat. 
Secondly,  although  it  is  possible  to  mix  graphics  and  text 
on  the  microcomputer  system,  it  is  not  done  in  a  straight  for- 
ward manner  and  utilizes  extra  memory.   For  memory  considera- 
tions and  convenience  it  was  decided  to  provide  sufficient 
information  for  interpretation  of  the  graphical  output 
below  each  graph.   It  is  felt  that  this  method  is  satisfac- 
tory and  creates  little  inconvenience  to  the  user.   To 
enhance  the  interpretation  of  the  graphical  output  much 
time  and  effort  was  devoted  to  developing  plotting  routines 
to  display  the  data  in  relation  to  axes  labeled  with  tic-marks 
and  boundaries  of  integer  vice  fractional  values. 

A  major  obstacle  that  was  overcome  was  the  identification 
and  correction  of  a  memory  management  problem  unique  to  the 
microcomputer  system  used.   This  problem  affected  only  those 
programs  requiring  graphical  output.   The  essence  of  the 
problem  was  due  to  the  size  of  the  programs  and  the  number 
of  variables  used,  parts  of  the  program  and  stored  variables 
were  being  over  written  in  memory  by  the  graphics  routines. 
This  problem  was  finally  solved  by  making  appropriate  changes 
in  the  memory  management  scheme.   For  a  more  detailed  explana- 
tion of  the  solution  see  Appendix  A. 
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III.   STATE  VARIABLE  PROGRAMS 

The  three  state  variable  programs  discussed  in  this  sec- 
tion are  modified  versions  of  the  programs  of  the  same  names 
found  in  [Ref.  1],   These  programs  may  be  used  as  tools  for 
the  analysis  and  design  of  linear  control  systems  expressed 
in  the  following  state  variable  form: 

x(t)  =  Ax(t)  +  bu(t) 
u(t)  =  K[r(t)  -  kTx(t)] 
y(t)  =  cTx(t) 

The  Basic  Matrix  program  (BASMAT)  described  in  part  A 
of  this  section  is  used  to  compute  the  determinant,  inverse, 
characteristic  polynomial,  and  eigenvalues  of  the  square 
matrix  A.   It  is  also  used  to  determine  the  resolvent  matrix 
(sI-A)~   and  state  transition  matrix  exp(At). 

The  Rational  Time  Response  program  (RTRESP),  described 
in  part  B,  is  used  to  determine  a  closed  form  expression 
for  the  time  response  of  a  system.   The  input  function  r(t) 
is  required  to  have  a  rational  time  response  and  no  repeated 
eigenvalues  are  allowed  in  the  combination  of  the  system 
and  input . 

The  Graphical  Time  Response  program  (RTRESP) ,  described 

in  part  C,  is  used  to  produce  a  graphical  display  of  the  time 

response  of  a  system  to  an  arbitrary  input. 
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RTRESP  and  GTRESP  can  be  used  to  study  open  loop  systems 
by  letting  K  equal  zero  and  unforced  systems  by  letting  r(t) 
equal  zero. 
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A.   BASMAT 

1.   Basic  Matrix  Program  (BASMAT) 

a.  Introduction 

When  done  by  hand,  matrix  manipulations  can  be 
quite  tedious  and  the  chances  of  an  error  being  made  are 
great.   In  the  study  of  linear  state  variable  analysis,  a 
computer  program  to  do  these  manipulations  is  almost 
essential.   It  can  do  the  necessary  manipulations  much  more 
quickly  and  accurately  and  allow  the  user  to  devote  his 
time  and  energy  to  design  and  analysis. 

b.  Description  of  Program 

BASMAT  [Ref.  1:   pp.  7,8]  takes  a  matrix  A, 
and  computes  the  determinant  of  A  (det  A),  the  inverse  of 
A  (A~  ),  the  characteristic  polynomial  (det  (sI_-A) )  ,  and 
eigenvalues  of  A  (Xt')as  well  as  the  resolvent  matrix 

<f>(s)    =  (sI-A)"1 
and  the  state  transition  matrix 

<j6(t)  =  exp(At) 
The  state  transition  matrix  is  expressed  as  matrix  coeffi- 
cients times  the  natural  modes  exp(\t)  with  the  complex 
eigenvalues  expressed  as  damped  sine  and  cosine  terms.   The 
resolvent  matrix  is  written  as 

(73(s)  =  adj(sI-A)/det(sI-A) 
and  the  matrix  numerator,  adj(sI-A),  is  output  as  matrix 
coefficients  of  powers  of  s  so  that  it  takes  the  form 
adj(sI-A)  =  F1  +  F2s  +  ...  +  F^-1 
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The  BASMAT  program  interactively  accepts  input 
of  a  matrix  and  calls  subroutines  to  perform  the  appro- 
priate calculations.   The  subroutines  used  are  CHREQ ,  CHREQA, 
PROOT,  DET,  SIMEQ,  and  STMST.   These  subroutines  are  listed 
below  with  a  brief  description. 

CHREQ.   This  subroutine  is  used  to  determine  the 
characteristic  polynomial  det(sI-A),  and  the  resolvent 
matrix  (sI.-A)-   for  the  matrix  A.   The  Leverrier  algorithm 

is  used  to  compute  the  resolvent  matrix 

-1 


0(s)  =  (sI-A) 


in  the  form  ^ 


N-i 


'     R.s 

0(s)  ■  — dT^T 


where  D(s)  is  the  characteristic  polynomial  det(sI-A). 

The  coefficients  of  the  characteristic  polynomial 
are  determined  by  subroutine  CHREQA. 

CHREQA.   This  subroutine  is  called  by  CHREQ  to 
determine  the  characteristic  polynomial,  det(sI-A),  for  the 
A  matrix.   The  principal-minor  method  is  used. 

PROOT.   This  subroutine  uses  a  modified  Bair- 
stow  method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.   It  is  used  to  calculate 
the  eigenvalues  of  the  A  matrix. 

DET.   This  subroutine  computes  the  determinant 
of  a  matrix.   A  gauss  elimination  method  to  place  the  matrix 
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in  upper  triangular  form  is  used.   It  is  called  by  CHREQA  to 
calculate  sub-determinants. 

SIMEQ.   This  subroutine  is  used  in  the  FORTRAN 
version  of  BASMAT  in  [Ref.  1:   pp.  7,8]  to  determine  the 
inverse  of  the  A  matrix.   In  the  BASIC  version  a  machine 
language  subroutine  [Ref.  2]  is  substituted  for  speed 
and  convenience. 

STMST.   This  subroutine  is  used  to  compute  the 
state  transition  matrix 

d>(t)    =  exp(At) 
for  a  matrix  A.   It  uses  the  Sylvester  Expansion  Theorem. 
Eigenvalues  of  the  A  matrix  must  be  provided.   This  routine 
can  not  handle  duplicate  eigenvalues. 

c.   Program  Translation  Problems 

See  section  II,  Translation  Considerations. 
2.   BASMAT  User's  Guide 

This  program  occupies  about  7k  bytes  of  memory  and 
does  not  utilize  the  graphics  pages,  therefore  it  is  not 
necessary  to  relocate  the  disk  operating  system. 

Push  "return"  after  each  input. 
STEP  1: 

SCREEN  PROMPT: 

BASIC  MATRIX  PROGRAM 
PROBLEM  IDENTIFICATION  -> 
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REMARKS : 

Label  the  problem.   Type  in  any  appropriate  combina- 
tion of  letters,  numbers,  and/or  symbols,  excluding 
commas.   This  input  is  limited  to  255  characters. 
STEP  2: 

SCREEN  PROMPT: 

PLANT  ORDER?  -> 
REMARKS : 

Enter  the  number  that  corresponds  to  the  order  of 
the  A  matrix.   Maximum  order  allowed  is  ten. 
STEP  3: 

SCREEN  PROMPT: 

INPUT  PLANT  MATRIX. 
A (row, column )= 
REMARKS : 

Input  the  A  matrix  by  typing  in  each  element  as 
prompted.   The  program  will  ask  for  each  element 
beginning  with  the  first  row  and  going  from  left  to 
right . 
STEP  4: 

SCREEN  PROMPT: 

HARDCOPY?  (Y/N)  -> 
REMARKS : 

After  this  step  is  completed  the  program  will  output 
the  A  matrix  for  reference.   This  will  be  followed  by 
the  inverse  of  the  A  matrix,  the  determinant  of  the  A 
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matrix,  the  matrix  coefficients  of  the  resolvent 
matrix  numerator,  the  characteristic  polynomial 
coefficients,  the  eigenvalues  of  the  plant  matrix, 
and  finally  the  elements  of  the  state  transition 
matrix.   See  Figure  1  for  a  sample  output. 
3.   BASMAT  Example 

This  example  will  use  the  system  matrix 


A  = 


0    10 

0    0    1 

-2   -3   -3 


to  demonstrate  the  use  of  the  BASMAT  program.   It  will  refer 

to  the  steps  described  in  the  previous  section,  BASMAT 

User's  Guide. 

Step  1.   Enter  "EXAMPLE  <CR>"  where  <CR>  =  return. 

Step  2.   Enter  "3  <CR>" 

Step  3.   Enter  "0  <CR>  1  <CR>  0  <CR>  0  <CR>  0  <CR>  1  <CR>  -2 

<CR>  -3  <CR>  -3  <CR>" 

Step  4.   Enter  "Y  <CR>" 

The  resulting  output  is  shown  in  Figure  1.   These 
results  are  interpreted  as  follows: 

The  plant  matrix  is  shown  just  as  it  was  entered  with 
rows  horizontal  and  columns  vertical.   This  is  the  format 
used  for  all  matrix  output. 

The  inverse  of  the  plant  matrix  is  shown  followed  by 
the  scalar  determinant  of  the  A  matrix. 
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The  matrix  coefficients  of  the  resolvent  matrix 
numerator  are  given  as  powers  of  s.  The  characteristic 
polynomial  is  listed  as  coefficients  of  powers  of  s. 

From  this  output ,  the  resolvent  matrix 
(j6(s)  =  adj(sI-A)/det(sI-A) 

may  be  written.   The  characteristic  polynomial,  det(sI-A), 

2   3 
is  2+3s+3s  +s  .   The  first  element  of  the  resolvent  matrix 

2 
numerator  is  3+3s+s   making  the  first  element  of  the 

resolvent  matrix 

2 


6(s)  = 


3+3s+s 


2+3s+3s2+s3 


The  real  and  imaginary  parts  of  the  eigenvalues  of 
the  plant  matrix  are  listed.   Finally,  the  elements  of  the 
state  transition  matrix  are  given.   The  first  element  of 
the  state  transition  matrix  can  be  written 

(f)fl(t)    =  0.333e"2t+0.667e"0,5t  cos  0  .  866t  +  l .  15e"°  ,5t  sin  0.866t 
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BASIC  MATRIX  PROGRAM 

PROBLEM  I  DENT I FI CAT  I  ON -EXAMPLE 

THE  PLANT  MATRIX 

8  1  9 

9  9  1 
-2            -3             -3 

THE  INVERSE  OF  THE  PLANT  MATRIX 

-1.5  -1.5  -.5 

1  9  9 

-9  1  9 

THE  DETERMINANT  OF  THE  PLANT  MATRIX 


THE  MATRIX  COEFFICIENTS  OF  THE  RESOLVENT  MATRIX  NUMERATOR 

THE  MATRIX  COEFFICIENT  OF  3A2 

1  9  9 

9  1  9 

9  9  1 

THE  MATRIX  COEFFICIENT  OF  SA 1 

3  1  9 

9  3  1 

-2  -3  9 

THE  MATRIX  COEFFICIENT  OF  Sv9 

3  3  1 

-2  9  0 

0-2  9 

********XX****£*****XXXXK*S£**£;(S  Figure  1   BASMAT  Output 
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Figure  1  (Cont . ) 

THE  CHARACTERISTIC  POLYNOMIAL  IN  ASCENDING  POWERS 

2   3   3   1 

THE  EIGENVALUES  OF  THE  PLANT  MATRIX 

REAL  IMAGINARY 

-2  9 

-.5  -.36602548  4 

-.5  .366025484 

******** ***X*****X***SX**** ****** 

THE  ELEMENTS  OF  THE  STATE  TRANSITION  MATRIX 


THE  MATRIX  COEFFICIENT  OF  EXP<-2)T 

.333333333    .333333333  .333333333 

- . 666666666       - . 666666666  - . 666666666 

1  .  33333333    1  . 33333333  1  . 33333333 

THE  MATR I X  COEFF I C I EMT  OF  EXP < - . 5) T* COS (  . 3668 25404) T 

.666666666         -.333333333  -.333333333 

.006666606     1  • 0000000/  .6  o 0006  O O O 

THE  MATR I X  COEFF I C I ENT  OF  EXP < - . 5) T*  S I N (  . 3668 2548 4) T 

1.15478054    1.73205031  .577358263 

-1.15478854   -.577358263  8 

3  -1.15478054  -.577350263 

********************************* 


Figure  1 
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B.   RTRESP 

1.   Rational  Time  Response  Program  (RTRESP) 

a.  Introduction 

Frequently  it  is  desirable  to  know  the  response 
of  a  system  as  a  function  of  time.   A  computer  program  can 
determine  this  quicker  and  more  accurately  than  by  hand. 

b.  Description  of  Program 

RTRESP  determines  the  time  response  in  closed 
form  of  the  closed  loop  system 

x(t )=Ax(t )+bu(t ) 

u(t)=K[r(t)-kTx(t)] 

y(t)=cTx(t) 

due  to  any  initial  conditions  x(0)  and  input  r(t)  for  t>0. 
The  system  must  have  a  rational  Laplace  transform  R(s)  with 
a  pole-zero  excess  of  at  least  one.   [Ref.  1:   pp.  11,12] 

The  input  r(t)  is  treated  by  forming  a  m   -order 
dynamic  system  whose  initial  condition  response  is  equal  to 
r(t)  for  a  specific  set  of  initial  conditions.   This  system 
is  combined  with  the  original  system  and  then  the  complete 
response  in  closed  form  is  determined  from  the  subroutine 
STMST.   The  order  of  the  combined  system  must  be  ten  or 
less. 

Various  primary  and  utility  subroutines  are  used 

in  RTRESP.   The  primary  subroutines  are  listed  below  with  a 

brief  description: 
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CHREQA.   This  subroutine  determines  the  charac- 
teristic polynomial,  det(sI-A),  for  the  matrix  A  using  the 
principal-minor  method. 

PROOT.  This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a  poly- 
nomial with  real  coefficients.  It  is  used  to  determine  the 
roots  of  the  numerator  and  denominator  of  R(s)  when  they  are 
entered  in  polynomial  form.  PROOT  is  also  used  to  determine 
the  eigenvalues  of  the  combined  system  matrix  (i.e.  the  roots 
of  the  characteristic  polynomial  determined  by  CHREQA). 

SEMBL.   This  subroutine  determines  the  coeffi- 
cients of  a  polynomial  from  the  roots  of  the  polynomial. 
It  is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  R(s)  when  they  are  entered 
in  the  factored  form.   This  subroutine  together  with  PROOT 
provide  the  feature  that  R(s)  may  be  entered  in  either 
factored  or  polynomial  form.   It  will  appear  in  the  output 
in  both  forms. 

STMST.   This  subroutine  computes  the  state  transi- 
tion matrix 

<j6(t)  =  exp(At) 
for  a  matrix  A.   It  uses  the  Sylvester  Expansion  Theorem. 
Eigenvalues  of  the  A  matrix  must  be  provided.   This  routine 
can  not  handle  duplicate  eigenvalues,  therefore  it  is 
necessary  that  the  combined  system  and  input  have  no 
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repeated  eigenvalues.   STMST  is  used  to  determine  the  state 

transition  matrix  of  the  augmented  system. 

c.   Program  Translation  Problems 

See  section  II,  Translation  Considerations. 

2.   RTRESP  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

RATIONAL  TIME  RESPONSE  PROGRAM  (RTRESP) 

THIS  PROGRAM  DETERMINES  THE  TIME  RESPONSE  OF  A 
CLOSED-LOOP  SYSTEM  DUE  TO  SPECIFIED  INITIAL 
CONDITIONS  AND  INPUT.   SYSTEM  MUST  HAVE  A  RA- 
TIONAL LAPLACE  TRANSFORM  WITH  A  POLE-ZERO  EXCESS 
OF  AT  LEAST  ONE. 

PROBLEM  IDENTIFICATION  -> 

REMARKS : 

Label  the  problem.   Type  in  any  appropriate  combina- 
tion of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.   This  input  is  limited  to  255 
characters.   This  is  also  the  name  used  to  save  and 
retrieve  the  disk  file  containing  the  problem 
description.   (See  steps  2  and  3) 
STEP  2: 

SCREEN  PROMPT: 

WAS  THIS  PROBLEM  DESCRIPTION  PREVIOUSLY 
SAVED?  (Y/N)  -> 

REMARKS : 

This  step  provides  the  option  of  recalling  a  pre- 
viously saved  problem  description  from  the  disk. 
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If  "Y"  is  typed,  the  program  looks  for  the  problem 

description  saved  under  the  label  typed  in  STEP  1. 

The  program  then  runs  to  completion  using  the 

retrieved  program  description.   If  an  "N"  is 

typed,  the  program  continues  with  STEP  3. 

STEP  3: 

SCREEN  PROMPT: 

DO  YOU  WANT  TO  SAVE  THIS  PROBLEM  DESCRIPTION? 
(Y/N)  -> 

REMARKS : 

A  positive  response  will  save  the  problem  descrip- 
tion, that  will  be  entered  in  the  following  steps, 
to  the  disk  under  the  label  entered  in  STEP  1. 
This  problem  may  be  retrieved  during  a  later  session 
by  entering  the  proper  label  in  STEP  1  and  typing 
"Y"  in  STEP  2. 
STEP  4: 

SCREEN  PROMPT: 

ORDER  OF  THE  SYSTEM?  -> 
REMARKS : 

The  total  order  of  the  system  and  the  input  must 
not  exceed  ten. 
STEP  5: 

SCREEN  PROMPT: 

INPUT  SYSTEM  (A)  MATRIX 
A(l,l)  = 
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REMARKS : 

Enter  the  elements  of  the  A  matrix  as  prompted. 

The  format  is  A( row, column )  =  . 

STEP  6: 

SCREEN  PROMPT: 

THE  A  MATRIX 

(display  of  the  A  matrix) 

ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

If  an  "N"  is  typed,  the  program  proceeds  to  STEP  9. 

To  correct  the  matrix  type  "Y"  and  the  program  will 

proceed  to  STEP  7. 

STEP  7: 

SCREEN  PROMPT: 

TYPE  ROW, COLUMN  OF  THE  ELEMENT  TO  BE 
CORRECTED  -> 

STEP  8: 

SCREEN  PROMPT: 

A  (row, column)  = 

REMARKS: 

Enter  the  correct  value.   The  program  will  return 

to  STEP  6. 

STEP  9: 

SCREEN  PROMPT: 

INPUT  THE  CONTROL  (B)  VECTOR 

B(l)  = 
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REMARKS : 

After  the  B  vector  is  entered,  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 
STEP  10: 

SCREEN  PROMPT: 

INPUT  THE  OUTPUT  (C)  VECTOR 
C(l)  = 
REMARKS : 

After  the  C  vector  is  entered,  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 
STEP  11: 

SCREEN  PROMPT: 

INPUT  THE  FEEDBACK  COEFFICIENTS 
FEEDBACK  COEFFICIENT  (1)  = 
REMARKS : 

After  the  feedback  coefficients  are  entered,  a 
correction  sequence  similar  to  steps  6  through  8 
is  encountered. 
STEP  12: 

SCREEN  PROMPT: 

INPUT  THE  GAIN  -> 
REMARKS : 

Enter  the  controller  gain. 
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STEP  13: 

SCREEN  PROMPT: 

INPUT  THE  INITIAL  CONDITIONS 

X0(1)  = 

REMARKS : 

After  the  initial  conditions  are  entered,  a 

correction  sequence  similar  to  steps  6  through 

8  is  encountered. 

STEP  14: 

SCREEN  PROMPT: 

ENTER  THE  INPUT  GAIN  -> 

REMARKS : 

Enter  the  gain  of  the  input  function,  R(s). 

STEP  15: 

SCREEN  PROMPT: 

■KEY1  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  NUMERATOR  POLYNOMIAL  -> 

REMARKS : 

Choose  the  preferred  method  of  entering  the  numerator 

polynomial  of  the  input  function,  R(s). 

STEP  16A: 

SCREEN  PROMPT  (assumes  coefficient  form  chosen): 

INPUT  POLYNOMIAL  COEFFICIENTS 

INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  S 

INPUT  COEFFICIENT  OF  S~0  -> 
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REMARKS : 

Enter  the  coefficient  along  with  its  algebraic 

sign  as  prompted.   The  program  assumes  that  the 

coefficient  of  the  highest  power  of  s  is  one. 

STEP  16B  (assumes  factored  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 

FOR  FACTOR  1 
REAL  PART  -> 

IMAGINARY  PART  -> 

REMARKS : 

Enter  the  factor  with  its  algebraic  sign.   Do  not 
enter  the  associated  root.  (e.g.  For  "(s-1)"  enter 
"-1".)  After  the  real  part  of  the  factor  is  entered, 
the  program  asks  for  the  imaginary  part,  allowing 
for  input  of  quadradic  factors.   If  a  non-zero 
imaginary  part  is  entered,  the  program  automatically 
enters  its  complex  conjugate. 
STEP  17: 

SCREEN  PROMPT: 

i KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR  POLYNOMIAL  -> 

REMARKS : 

The  denominator  is  entered  just  as  described  for 

the  numerator  in  steps  15  and  16. 
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After  this  step  is  completed,  the  program  will 
output  all  input  for  reference  followed  by  the 
time  response  of  x(t)  and  the  system  output, 
y(t). 

3.   RTRESP  Example 

It  is  desired  to  know  the  time  response  of  the 
closed  loop  system 


x(t)  = 


0    1 
-1   -1 


x(t)   + 


u(t) 


u(t)  =  5{r(t)  -  [1.0   0.5]  x(t)} 


y(t)  =  [1   1]  x(t) 


if  the  input  function  is  given  by 

R(s)  =  2.0  (1/s) 
and  the  initial  conditions  are 

x(t)  =  0 
Referring  to  the  steps  described  in  the  pervious  sec- 
tion, RTRESP  User's  Guide,  this  problem  would  be  entered  as 
follows : 

Step  1.   Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 
return . 

Step  2  and  Step  3.   Enter  "N  <CR> . " 
Step  4.   Enter  "2  <CR>." 

Step  5.   Enter  "0  <CR>  1  <CR>  -1  <CR>  -1  <CR>." 
Step  6.   Enter  "N  <CR> . " 
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Step  9.   Enter  "0  <CR>  1  <CR> . "   Enter  appropriate 
response  for  correction  sequence. 

Step  10.   Enter  "1  <CR>  1  <CR> . "   Enter  appropriate 
response  for  correction  sequence. 

Step  11.   Enter  "1  <CR>  .5  <CR>."   Enter  appropriate 
response  for  correction  sequence. 

Step  12.   Enter  "5  <CR>." 

Step  13.   Enter  "0  <CR>  0  <CR> . "   Enter  appropriate 
response  for  correction  sequence. 

Step  14.   Enter  "2  <CR>." 

Step  15.   Enter  "P,0  <CR>." 

Step  16A.   Enter  "1  <CR>." 

Step  17.   Enter  *'F,1  <CR>  0  <CR>." 
The  program  will  now  run  to  completion  resulting  in  the  out- 
put shown  in  Figure  2. 

The  output  begins  with  the  program  name  and  problem 
identification.   The  A,  B,  and  £  matrices  are  shown  for 
reference  followed  by  the  feedback  coefficients  and  the 
controller  gain.   Next,  the  initial  conditions  and  input 
gain  are  listed.   The  numerator  and  denominator  polynomials 
of  the  input  function  are  given  in  both  polynomial  coeffi- 
cient form  and  factored  form  regardless  of  the  method  in 
which  they  were  entered.   The  polynomial  coefficient  form  is 
given  as  a  list  of  coefficients  from  left  to  right  in 
ascending  powers  of  s.   The  highest  power  of  s  is  always 
one.   The  polynomial  factored  form  appears  as  a  list  of  the 
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real  and  imaginary  parts  of  each  root  of  the  polynomial. 
Notice  that  the  numerator  is  only  given  in  coefficient 
form  since  it  is  of  zero  order  and  no  roots  exist . 

The  time  response  of  the  state  x(t)  is  given  in  the 
form  of  vector  coefficients  of  the  various  natural  modes  of 
the  system  and  the  input  function.   From  this  output 
x(t)  is  seen  to  be 

x1(t)  =  -1.67e  x''OL   cos  1.71t  -  1 . 70e  x# 'ox   sin  1.71t  +  1.67 

-1  75t 
x2(t)  =  5.83e  x*'ot   sin  1.71t 

The  time  response  of  the  output  y(t)  is  given  as  scalor 
coefficients  of  the  same  natural  modes  as  x(t)  so  it  is  seen 
that 

-1  7c»t  -1  7cit 

y(t)  =  -1.67e  X''DZ      cos  1.71t  +  4.13e  x-'OL   sin  1.71t  +  1.67 
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RATIONAL  TIME  RESPONSE 

PROBLEM  IDENTIFICATION  ->  EXAMPLE 

THE  A  MATRIX 

8  1 

-1     -1 

THE  B  MATRIX 

9  1 

THE  C  MATRIX 

1      1 

FEEDBACK  COEFFICIENTS 

1      .5 

GAIN  =  5 

INITIAL  CONDITIONS,  X<9> 
9      9 
RGAIN  =  2 


NUMERATOR  POLYNOMIAL  OF  R<S)  -  ASCENDING  PWRS  OF  S 

1 


Figure  2   RTRESP  Output 
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Figure  2  (Cont . ) 

DENOMINATOR  POLYNOMIAL  OF  R<S)  -  ASCENDING  PWRS  OF  3 

8      1 

DENOMINATOR  ROOTS  ARE 

REAL  PART        IMAGINARY  PART 

0  0 

THE  TIME  RESPONSE  OF  THE  STATE  X(T) 

THE  VECTOR  COEFF I C I ENT  OF  EXP  < -  1 . 75) T *  COS  <  1  . 7 1 39 1365) T 
-1.6666666?  9 

THE  VECTOR  COEFF I C I ENT  OF  EXP < -  1 . 75) T* S I N< 1 . 7 1 39 1 365) T 
-1.78  175324      5.33459966 

THE  VECTOR  COEFFICIENT  OF  EXP03)T 
1.66466667       ? . 3  1322575E-  19 

THE  TIME  RESPONSE  OF  THE  OUTPUT  Y(T) 

THE  COEFFICIENT  OF  EXP< -  1  . 75) T*COS<  1  . 7  13?  1365) T 

-  1  . 66666667 

THE  COEFF I C I ENT  OF  EXP ( - 1 . 75) T* S I N <  1  . 7  1 39  1 365) T 

4.  13234143 

THE  COEFFICIENT  OF  EXP03)T 

1  . 66666667 


Figure  2 
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C.   GTRESP 

1.   Graphical  Time  Response  Program  (GTRESP) 

a.  Introduction 

Knowing  the  response  of  a  system  as  a  function  of 
time  can  be  very  helpful  in  the  study  and  analysis  of  that 
system.   The  graphical  display  of  this  response  can  give 
much  insight  into  a  system's  response  characteristics. 

b.  Description  of  Program 

GTRESP  [Ref.  1:   pp.  22-28]  determines  and 
graphically  displays  the  time  response  of  the  closed  loop 
system 

x(t)  =  Ax(t)  +  bu(t) 
u(t)  =  K  [r(t)  -  kT  x(t)] 
y(t)  =  cT  x(t) 

due  to  any  initial  conditions  x(0)  and  input  r(t). 

The  basic  purpose  of  this  program  is  similar 
to  that  of  the  RTRESP  program  described  earlier.   The 
difference  between  the  two  programs  is  that  GTRESP  will 
determine  the  time  response  for  arbitrary  input  functions 
which  may  not  have  rational  Laplace  transforms  but  RTRESP 
requires  a  rational  Laplace  transform.   Also  GTRESP  produces 
a  graphical  display  of  the  time  response  instead  of  a 
closed  form  expression  as  in  the  RTRESP  program. 

The  GTRESP  program  uses  a  fourth-order  Runge-Kutta 
numerical  integration  algorithm  to  calculate  the  time  response 
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The  main  subroutines  used  are  CALCU,  RUNGE ,  TRESP ,  and  YDOT. 
Also  various  utility  and  plotting  subroutines  are  used. 

CALCU.   This  subroutine  is  called  by  TRESP  to 
determine  the  reference  input  r(t)  and  the  control  input 

u(t)=K  [r(t)  -  kTx(t)] 

The  reference  input  r(t)  must  be  defined  by  the  user  by 
inserting  the  appropriate  BASIC  coding  into  the  subroutine 
between  line  numbers  5010  and  5500.   The  reference  input  is 
represented  by  the  variable  R  and  the  control  input  by  the 
variable  U. 

RUNGE.   This  subroutine  is  called  by  TRESP  and 
contains  the  actual  fourth-order  Runge-Kutta  integration 
algorithm.   It  must  be  executed  four  times  for  each  inte- 
gration step. 

TRESP.   This  subroutine  is  the  driving  subroutine 
which  calls  the  subroutines  CALCU,  RUNGE,  and  YDOT  along  with 
the  necessary  plotting  routines.   It  calculates  the  time 
response  of  the  closed  loop  linear  system  described  by  the 
input  parameters  and  plots  the  desired  variables. 

YDOT.   This  subroutine  is  called  by  TRESP  to 
compute  the  derivative 
x(t)=Ax(t)+bu(t) 

It  is  designed  to  handle  linear  systems  but  can  easily  be 
modified  to  handle  nonlinear  and  time  varying  systems  which 
would  give  GTRESP  a  nonlinear  and  time  varying  capability. 
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c.   Program  Translation  Problems 

See  section  II,  Translation  Considerations. 

2.   GTRESP  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

GRAPHICAL  TIME  RESPONSE  PROGRAM  (GTRESP) 

THIS  PROGRAM  DETERMINES  AND  GRAPHICALLY 
DISPLAYS  THE  TIME  RESPONSE  OF  A  CLOSED- 
LOOP  SYSTEM  DUE  TO  SPECIFIED  INITIAL 
CONDITIONS  AND  INPUT. 

PROBLEM  IDENTIFICATION  -> 

REMARKS : 

Label  the  problem.   Type  in  any  appropriate 
combination  of  letters,  numbers,  and/or  symbols, 
excluding  commas  and  colons.   This  input  is 
limited  to  255  characters.   This  is  also  the  name 
used  to  save  and  retrieve  the  disk  file  containing 
the  problem  description.   (See  steps  2  and  3) 
STEP  2: 

SCREEN  PROMPT: 

WAS  THIS  PROBLEM  DESCRIPTION  PREVIOUSLY 
SAVED?  (Y/N)  -> 

REMARKS : 

This  step  provides  the  option  of  recalling  a  pre- 
viously saved  problem  description  from  the  disk. 
If  "Y"  is  typed  the  program  looks  for  the  problem 
description  saved  under  the  label  typed  in  STEP  1. 
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The  program  then  runs  to  completion  using  the 

retrieved  program  description.   If  an  "N"  is 

typed  the  program  continues  with  STEP  3. 

STEP  3: 

SCREEN  PROMPT: 

DO  YOU  WANT  TO  SAVE  THIS  PROBLEM  DESCRIPTION? 
(Y/N)  -> 

REMARKS : 

The  problem  description  will  be  entered  in  the 
following  steps. 

A  positive  response  will  save  the  problem  des- 
cription to  the  disk  under  the  label  entered  in 
STEP  1.   This  problem  may  be  retrieved  during  a 
later  session  by  entering  the  proper  label  in  STEP 
1  and  typing  "Y"  in  STEP  2. 
STEP  4: 

SCREEN  PROMPT: 

ORDER  OF  THE  SYSTEM?  -> 
REMARKS : 

The  total  order  of  the  system  and  the  input  must  not 
exceed  ten. 
STEP  5: 

SCREEN  PROMPT: 

INPUT  SYSTEM  (A)  MATRIX 
A(l,l)  = 
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REMARKS : 

Enter  the  appropriate  element  of  the  A  matrix.   The 

format  is  A(row, column )  =  . 

STEP  6: 

SCREEN  PROMPT: 

THE  A  MATRIX 

(display  of  the  A  matrix) 

ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

If  an  "N"  is  typed,  the  program  proceeds  to  STEP  9. 

To  correct  the  matrix  type  "Y"  and  the  program  will 

proceed  to  STEP  7. 

STEP  7: 

SCREEN  PROMPT: 

TYPE  ROW, COLUMN  OF  THE  ELEMENT  TO  BE 
CORRECTED  -> 

STEP  8: 

SCREEN  PROMPT: 

A (row, column )  = 

REMARKS : 

Enter  the  correct  value.   The  program  will  return  to 

STEP  6. 

STEP  9: 

SCREEN  PROMPT: 

INPUT  THE  CONTROL  (B)  VECTOR 

B(l)  = 
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REMARKS : 

After  the  B  vector  is  entered  a  correction 
sequence  similar  to  steps  6  through  8  is 
encountered. 
STEP  10: 

SCREEN  PROMPT: 

INPUT  THE  OUTPUT  (C)  VECTOR 
C(l)  = 
REMARKS : 

After  the  C  vector  is  entered  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 
STEP  11: 

SCREEN  PROMPT: 

INPUT  THE  FEEDBACK  COEFFICIENTS 
FEEDBACK  COEFFICIENT  (1)  = 
REMARKS : 

After  the  feedback  coefficients  are  entered  a  correc- 
tion sequence  similar  to  steps  6  through  8  is 
encountered. 
STEP  12: 

SCREEN  PROMPT: 

INPUT  THE  CONTROLLER  GAIN  -> 
REMARKS : 

Enter  the  controller  gain. 
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STEP  13: 

SCREEN  PROMPT: 

INPUT  THE  INITIAL  CONDITIONS 
X0(1)  = 
REMARKS : 

After  the  initial  conditions  are  entered  a  correc- 
tion sequence  similar  to  steps  6  through  8  is 
encountered. 
STEP  14: 

SCREEN  PROMPT: 

ENTER  THE  FOLLOWING  PARAMETERS: 
NOTE. . (TF-TZ)/(DT*FR)  =<  100 

INITIAL  TIME  (TZ)  -> 

REMARKS: 

Enter  the  initial  time  of  the  time  interval  of 

interest.   Due  to  program  constraints  the  initial 

and  final  times  and  the  time  step  and  frequency  of 

output  must  be  chosen  so  that  they  satisfy  the 

relation . 

(TF-TZ)  /  (DT*FR)  =<  100 

STEP  15: 

SCREEN  PROMPT: 

FINAL  TIME  (TF)  -> 

REMARKS : 

Enter  the  final  time  of  the  time  interval  of 

interest . 
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STEP  16: 

SCREEN  PROMPT: 

TIME  STEP  (DT)  -> 
REMARKS : 

Enter  the  time  increment  for  each  step. 
STEP  17: 

SCREEN  PROMPT: 

FREQUENCY  OF  OUTPUT  (FR)  -> 
REMARKS : 

Enter  an  integer  n  and  the  program  will  print  out 
data  on  every  nth  time  step  iteration.   Ensure 
that  the  relation  expressed  in  the  remarks  of  STEP 
14  is  satisfied. 
STEP  18: 

SCREEN  PROMPT: 

YOU  MAY  PLOT  UP  TO  8  VARIABLES  VS  TIME. 
VARIABLE   NUMBER     VARIABLE   NUMBER 


XKT) 

1 

X8(T) 

8 

X2(T) 

2 

X9  ( T ) 

9 

X3(T) 

3 

X10(T) 

10 

X4(T) 

4 

E(T) 

11 

X5(T) 

5 

U(T) 

12 

X6(T) 

6 

Y(T) 

13 

X7(T) 

7 

R(T) 

14 

HOW  MANY  VARIABLES  TO  PLOT?  MAX=8  "> 
REMARKS : 

You  may  plot  up  to  8  variables  vs  time.   These 
variables  will  be  plotted  on  the  same  plot  in  the 
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output.   Enter  the  number  of  variables  you  want  to 

appear  on  the  plot. 

STEP  19: 

SCREEN  PROMPT: 

TYPE  THE  VARIABLE  NUMBER  <CR>  -> 

REMARKS : 

Type  the  number  associated  with  the  variable  of 

interest.   After  carriage  return  <CR>  is  typed,  the 

program  will  continue  to  prompt  for  another  variable 

number  until  the  number  of  variables  that  the  user 

indicated  in  STEP  18  has  been  entered. 

STEP  20: 

SCREEN  PROMPT: 

YOU  CHOSE  THE  FOLLOWING  VARIABLES: 

(list  of  variable  numbers  entered  in  STEP  19) 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

To  make  a  change  in  the  variables  to  be  plotted 

type  a  "Y"  and  the  program  will  return  to  STEP 

18.   If  an  "N"  is  typed,  the  program  will  run  to 

completion . 

3.   GTRESP  Example 

It  is  desired  to  determine  and  plot  the  time  response 

of  the  error 

e(t)=r(t)-y(t) 

the  input  r(t),  and  the  state  variable  x2(t)  for  the  system 
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x(t)  = 
u(t)  = 

y(t)  = 


o   l 
-l  -l 


x(t)  + 


u(t) 


5{r(t)  -  [1.0   0.5]  x(t)} 
[1   1]  x(t) 


The  time  interval  of  interest  is 

0<t<4 
The  iteration  step  size  (DT)  is  chosen  to  be  0.01  and  printed 
output  is  desired  every  ten  steps  so  FREQ  equals  10.   The 
initial  conditions  are 

x(t)=0 
and  the  input  function  is 

r(t)  =  5.0     if  0<t<l 

r(t)  =  0       otherwise. 

In  order  to  define  this  input  function  the  following 

BASIC  coding  is  inserted  between  line  numbers  5010  and  5500 

in  the  subroutine  CALC. 

5020  IF  TZ  >  1  GOTO  5040 

5030  R  =  5:   GOTO  5500 

5040  R  =  0 

5500  REM   END  OF  ROUTINE  DESCRIBING  R(T) 

Referring  to  the  steps  described  in  section  2, 
GTRESP  User's  Guide,  this  problem  would  be  entered  as 
follows: 

Step  1.   Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 
return. 

Step  2  and  Step  3.   Enter  "N  <CR> . " 

Step  4.   Enter  "2  <CR>." 
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Step  5.   Enter  "0  <CR>  1  <CR>  -1  <CR>  -1  <CR> . " 

Step  6.   Enter  "N  <CR>." 

Step  9.   Enter  "0  <CR>  1  <CR> . "   Enter  appropriate 
response  for  correction  sequence. 

Step  10.   Enter  "1  <CR>  1  <CR>."   Enter  appropriate 
response  for  correction  sequence. 

Step  11.   Enter  "1  <CR>  .5  <CR>."   Enter  appropriate 
response  for  correction  sequence. 

Step  12.   Enter  "5  <CR>." 

Step  13.   Enter  "0  <CR>  0  <CR>."   Enter  appropriate 
response  for  correction  sequence. 

Step  14.   Enter  "0  <CR>." 

Step  15.   Enter  "4  <CR> . " 

Step  16.   Enter  ".01  <CR> . " 

Step  17.   Enter  "10  <CR>."   Ensure  that  the  constraint 
given  in  step  14  of  GTRESP  User's  Guide  is  met.   That  con- 
straint is  (TF-TZ)/(DT*FR)  =   100 
In  this  case 
(4-0)/( .01*10)=40 
so  the  constraint  is  met. 

Step  18.   Enter  "3  <CR>." 

Step  19.   Enter  "11  <CR>  14  <CR>  2  <CR> . "   This  causes 
the  error  e(t),  the  input  r(t),  and  the  state  variable  x~(t), 
respectively,  to  be  plotted  as  desired  from  the  problem 
statement . 
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The  program  will  now  run  to  completion  resulting  in 
the  output  shown  in  Figure  3. 

The  output  begins  with  the  program  name  and  problem 
identification.   The  A,  B,  and  C  matrices  are  shown  for 
reference  followed  by  the  feedback  coefficients  and  the 
controller  gain.   Next  the  initial  conditions  and  time 
parameters  are  listed. 

The  second  page  of  output  lists  in  tabular  form  the 
value  of  time  t  and  the  corresponding  values  of  the  output 
y(t),  the  control  u(t),  and  all  of  the  state  variables. 
The  user  has  no  control  over  the  variables  output  in  this 
form. 

The  third  page  of  output  is  the  graph  itself.   Below 
the  graph  is  enough  information  to  properly  interpret  the 
plotted  data. 
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GRAPHICAL  TIME  RESPONSE 

PROBLEM  IDENTIFICATION  ->  EXAMPLE 

THE  A  MATRIX 


0 
-1 

1 
-1 

THE  B 

MATRIX 

9 

1 

THE  C 

MATRIX 

1 

1 

FEEDBACK  COEFFII 

:  I  ENT 

1 

.5 

GAIN  = 

=  5 

INITIAL  C ON D I T I  ON S ,  X ( 9 ) 

9  0 

INITIAL  TIME  =  0  FINAL  TIME  =  4 

TIME  STEP  =  .0  1  FREQUENCY  OF  OUTPUT  =  10 

.-%  .-\  -■>  A  A  A  .«\  j-s  ,-\  *\  ,-\  *\  A  s:  ••%  •■%  A  *\  v\  /•.  /«.  .■■.  A  .-*.  A   .-\  /s  A  .-v  A  A  .■*.  A  A  A  t:  •  %  ••%  *\  «-%  .■•,  A  A  •  %  s\  A  .'•.  «■%  .•-.  .•% 


Figure  3   GTRESP  Output 


47 


Figure    3    (Cont .  ) 

GRAPHICAL    TIME    RESPONSE 

PROBLEM    IDENTIFICATION    ->    EXAMPLE 


Y(T) 


IKT) 


XKT) 


:2<T) 


.1 

.2 
.3 

.4 

.5 

.599999998 

.699999996 

.799999994 

.899999991 

.999999939 

1.89999999 

1.19999998 

1.29999998 

1.39999998 

1.49999998 

1.59999998 

1.69999998 

1.79999997 

1.39999997 

1.99999997 

2.89999996 

2.19999996 

2.29999996 

2.39999996 

2.49999996 

2.59999995 

2.69999995 

2.79999995 

2.89999995 

2.99999995 

3.89999994 

3.19999994 

3.29999994 

3.39999994 

3.49999993 

3.59999993 

3.69999993 

3.79999993 

3.89999993 

3.99999992 


2.19944778 

3.34346675 

5.8265394 

5.31339897 

6.28493982 

6.51852439 

6.55142585 

6.46812831 

6.28025128 

6.84733255 

3.62837346 

1.69868535 

.258794762 

-.768353826 

-1.4543845 

-1.36487253 

-2.35893263 

-2.89253483 

-2.31880576 

-1.35192291 

-1.64667353 

-1.41988913 

-1.18722972 

-.964837842 

-.758194579 

-.574947468 

-.416853781 

-.284448821 

-.176816196 

-.8928623427 

-.8276993776 

.3198515208 

.3518366439 

.3718853287 

.381588273 

.3847905445 

.832838823 

.8772893216 

.8694859387 

.8684867363 


25 

19.2237187 

14.3947762 

18.4766595 

7.39775741 

5.3653327 

3.37643244 

2.22689887 

1.51321108 

1.14451516 

1.83693376 

-18.1932386 

-13.1375299 

-8.91382349 

-5.43739354 

-2.79648681 

-.759576686 

.734871188 

1.71656348 

2.33653976 

2.65573162 

2.7468357 

2.67275259 

2.48587785 

2.22396883 

1.93564372 

1.6314455 

1.33489664 

1.35867972 

.318784287 

.595121516 

.413224971 

.264229347 

. 145923933 

.8552849217 

-.8115845112 

-.85793434 

-.8877345728 

-.184315806 

-.118744482 

-.189696615 


.111864751 

.39362276 

.782796737 

1.22749837 

1.6389279 

2.13393263 

2.55813361 

2.93459526 

3.26194266 

3.53319394 

3.65693377 

3.5564866 

3.38641664 

2.96369324 

2.57239891 

2.16798317 

1.77333415 

1.48593864 

1.37618985 

.73963826 

.547919246 

.34993815 

.  192873584 

.3724495393 

-.3168649382 

-.8776387299 

-.117134954 

-.139822866 

-.147465518 

-.145986564 

-.137598611 

-.124743259 

-.189486217 

-.3938872894 

-.3768934685 

-.3616168885 

-.8477441939 

-.3355633192 

-.8251381781 

-.3166881488 


2.33333303 

3.454344 

4.24374261 

4.5859989 

4.59681112 

4.37162176 

3.99323644 

3.5i55iw85 

3.31838362 

2.58333861 

-.8365653122 

-1.35738124 

-3.34762137 

-3.73255737 

-4.82728341 

-4.33197569 

-3.33223678 

-3.49844268 

-3.38699561 

-2.64155317 

-2.19459277 

-1.76987733 

-1.33018331 

-1.33648655 

-.742129671 

-.497316738 

-.299748746 

-.145425754 

-.3293586731 

.353924521 

.189391233 

.  14379473 

.168442361 

.16489261 

.  158398742 

. 146497353 

.138582217 

.112352641 

.8946741163 

.3778949276 


Figure  3 
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Figure  3  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  GRAPHICAL  TIME  RESPONSE  XX 


ABSCISSA  ->  TIME  AXIS 
ORDINATE  ->  RESPONSE  MAGNITUDE 
TIC  MARKS  SHOW  INTERVALS  OF  UNITY 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  9  TO  4 

ORDINATE,  -5  TO  5 

SYSTEM  RESPONSE 

VARIABLE         SYMBOL 
X2  + 
ERRORX 
INPUTT 


Figure  3 
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IV.   TRANSFER  FUNCTION  PROGRAMS 

The  two  transfer  function  programs  discussed  in  this 
section  are  modified  versions  of  programs  of  the  same  names 
found  in  [Ref .  1] . 

The  Frequency  Response  program  (FRESP),  discussed  in 
part  A  of  this  section,  determines  and  plots  the  frequency 
response  of  a  transfer  function  over  a  specified  range  of 
frequencies.   The  output  may  take  the  form  of  rectangular 
Bode  plots  or  a  polar  Nyquist  plot  or  both  as  desired. 

The  Root  Locus  program  (RTLOC),  discussed  in  part  B, 
calculates  and  plots  the  root  locus  of  a  transfer  function 
for  a  specified  range  of  gains.   It  is  also  possible  to 
enlarge  a  small  rectangular  section  of  the  root  locus  for 
more  detail. 
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A.   FRESP 

1.   Frequency  Response  Program  (FRESP) 

a.  Introduction 

The  response  of  a  system  as  a  function  of  fre- 
quency is  a  very  important  characteristic  of  that  system. 
Some  common  ways  of  graphically  displaying  the  frequency 
response  of  a  system  include  the  use  of  amplitude  and  phase 
Bode  plots  and  Nyquist  diagrams.   These  are  valuable  tools 
for  system  analysis. 

b.  Description  of  Program 

The  FRESP  program  [Ref.  1:   pp.  105-113]  is 
used  to  determine  the  frequency  response  of  a  rational  trans- 
fer function  G(s)  of  the  form 

G(s)  =  K-^ 


B(s) 

where 

a/  n  M-l     M 

A(s)  =  a-  +  a0s  +  ...  +  a„s    +  s 
12  M 

N-l     N 
B(s)  =  b-  +  bns  +  ...  +  b„s     +  s 
12  N 

The  output  may  take  the  form  of  Bode  plots  or  a  Nyquist  dia- 
gram or  both  in  addition  to  tabular  data. 

FRESP  gives  the  user  the  option  of  supplying 
discrete  frequency  values  or  allowing  the  program  to  linearly 
or  logarithmically  interpolate  frequency  values  between  two 
limit  values.   The  complex  number  G(jCJ)  can  be  computed  for 
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each  frequency  value.   G(j(jJ)  can  be  written  in  rectangular 
form  as 

G(jGJ)  =  R(6l))  +  JX(GJ) 
where  R( CO)    and  X(QJ)    are  real  values.   The  magnitude  and 
phase  of  G(  j  (jj)    can  then  be  written  as 

G(JCU)  =  [R2(CJ)  +  X2(GJ)]* 
arg  G(JCJ)  =  arctan  X(OJ)/R(0J) 

The  FRESP  program  uses  the  subroutines  PROOT , 
PVAL,  and  SEMBL.   The  subroutines  MAXI ,  PHNOM,  GRAPH,  and 
SPLIT  that  appear  in  the  FORTRAN  version  do  not  appear  in  the 
BASIC  version.   They  are  incorporated  into  the  main  program 
and  into  various  plotting  subroutines  written  for  the  specific 
microcomputer  system  used. 

PROOT.   This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.   It  is  used  to  determine 
the  roots  of  the  numerator  and  denominator  of  G(s)  when 
they  are  entered  in  polynomial  form. 

SEMBL.   This  subroutine  determines  the  coeffi- 
cients of  a  polynomial  from  the  roots  of  the  polynomial. 
It  is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  G(s)  when  they  are  entered  in 
the  factored  form.   This  subroutine  together  with  PROOT 
provides  the  feature  that  G(s)  may  be  entered  in  either 
factored  or  coefficient  form  and  it  will  appear  in  the  output 

in  both  forms. 
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PVAL.   This  subroutine  is  used  to  evaluate  a 
polynomial  A(s)  with  real  coefficients  for 
s  =  PR  +  jPI 

PVAL  is  executed  twice  for  each  frequency  value  used,  once 
for  the  numerator  of  the  transfer  function  and  again  for 
the  denominator. 

c.   Program  Translation  Problems 

In  addition  to  the  programing  considerations  dis- 
cussed in  section  II,  an  output  anomaly  was  traced  to  an 
apparent  oversight  in  the  subroutine  PVAL.   It  was  discovered 
that  a  zero  order  numerator  over  a  second  order  denomina- 
tor with  a  free  s  such  as 

l/s(s+10) 
was  plotted  as 

l/(s+10) 

The  reason  for  this  was  that  PVAL  treated  a  zero  order 

polynomial  and  a  first  order  polynomial  both  as  a  first  order 

polynomial  which  in  this  case  resulted  in  the  free  s  in 

the  denominator  being  cancelled.   The  reason  for  this  can 

be  seen  by  examining  the  portion  of  FORTRAN  code  from  PVAL 

in  [Ref.  1:   p.  164]  repeated  below: 

P=CMPLX(A(NN+1) ,0. ) 
DO  100  J=1,NN 
100  P=P*S+A(NN+1-J) 
VR=REAL(P) 

In  this  portion  of  code  NN  represents  the  order 

of  the  polynomial  being  evaluated  and  the  intention  is  that 

the  DO  loop  be  executed  a  number  of  times  equal  to  the  order 
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of  the  polynomial.   But  because  of  the  nature  of  the  DO  loop, 

if  the  polynomial  is  of  zero  order  (i.e.  NN  equals  zero) 

the  loop  will  be  executed  once  just  as  if  the  polynomial 

was  of  first  order  [Ref.  3].   One  way  to  correct  this 

problem  is  to  modify  the  section  of  PVAL  as  shown  here: 

P=CMPLX(A(NN+1) ,0. ) 
IF  (NN)  110,110,90 
90  DO  100  J=1,NN 
100  P=P*S+A(NN+1-J) 
110  VR=REAL(P) 

The  essence  of  the  modification  is  to  skip  the  DO  loop  if 
the  polynomial  has  order  zero.   Translating  this  modifica- 
tion into  the  BASIC  version  of  FRESP  solved  the  problem. 
2.   FRESP  User's  Guide 

Push  "return"  after  each  input. 
STEP  1: 

SCREEN  PROMPT: 

FREQUENCY  RESPONSE  PROGRAM  (FRESP) 

THIS  PROGRAM  OBTAINS  AND  PLOTS  THE  FREQUENCY 
OF  A  RATIONAL  TRANSFER  FUNCTION  OVER  A  SPECIFIED 
RANGE  OF  FREQUENCIES.   BOTH  RECTANGULAR  BODE 
PLOTS  AS  WELL  AS  A  POLAR  NYQUIST  PLOT  CAN  BE 
OBTAINED. 

PROBLEM  IDENTIFICATION  -> 
REMARKS : 

Label  the  problem.   Type  in  any  appropriate  combina- 
tion of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.   This  input  is  limited  to  255 
characters . 
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STEP  2: 

SCREEN  PROMPT: 

INPUT  THE  GAIN  -> 

REMARKS : 

Enter  the  transfer  function  gain. 

STEP  3: 

SCREEN  PROMPT: 

•KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  NUMERATOR  ,A(S)  -> 
REMARKS : 

Choose  the  preferred  method  of  entering  the  numera- 
tor polynomial  of  the  transfer  function,  G(s). 
STEP  4A  (assumes  coefficient  form  chosen): 
SCREEN  PROMPT: 

INPUT  POLYNOMIAL  COEFFICIENTS 
INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  S 
INPUT  COEFFICIENT  OF  S~0  -> 
REMARKS : 

Enter  the  coefficient  along  with  its  algebraic  sign 
as  prompted.   The  program  assumes  that  the  coeffi- 
cient of  the  highest  power  of  s  is  one. 
STEP  4B  (assumes  factored  form  chosen): 
SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 

FOR  FACTOR  1 
REAL  PART  -> 
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IMAGINARY  PART  -> 

REMARKS : 

Enter  the  factor  with  its  algebraic  sign.   Do  not 
enter  the  associated  root.   (e.g.  For  "(s-1)" 
enter  "-1".)   After  the  real  part  of  the  factor  is 
entered,  the  program  asks  for  the  imaginary  part 
allowing  for  input  of  quadradic  factors.   If  a 
non-zero  imaginary  part  is  entered,  the  program 
automatically  enters  its  complex  conjugate. 
STEP  5: 

SCREEN  PROMPT: 

■KEY1  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
'KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR,  B(S)  -> 

REMARKS : 

The  denominator  is  entered  just  as  described  for 

the  numerator  in  steps  3  and  4. 

STEP  6: 

SCREEN  PROMPT: 

MINIMUM  FREQUENCY  VALUE  -> 

REMARKS : 

Enter  the  minimum  frequency  value  in  radians  per 

second  of  the  frequency  range  of  interest. 

STEP  7: 

SCREEN  PROMPT: 

MAXIMUM  FREQUENCY  VALUE  -> 
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REMARKS : 

Enter  the  maximum  frequency  value  in  radians  per 
second  of  the  frequency  range  of  interest. 
STEP  8: 

SCREEN  PROMPT: 

NUMBER  OF  FREQUENCY  VALUES  TO  BE  USED  -> 

REMARKS : 

Enter  an  integer  from  1  to  200.   Since  the  horizon- 
tal resolution  capability  of  the  graphics  on  the 
microcomputer  system  used  is  191  pixels,  any  number 
of  values  greater  than  191  adds  no  more  detail  to  the 
Bode  plots.   For  the  Nyquist  diagram  the  greater  the 
number  of  values  the  greater  the  detail  due  to  the 
nature  of  the  plot.   If  a  list  of  discrete  fre- 
quency values  will  be  supplied  (i.e.  option  1  in 
step  9)  enter  the  number  of  frequencies  in  that 
list . 
STEP  9: 

SCREEN  PROMPT: 

0  =  LOGARITHMIC  INTERPOLATION 

1  =  DISCRETE  VALUES  SUPPLIED 

2  =  LINEAR  INTERPOLATION 

CHOOSE  ONE  -> 
REMARKS : 

Make  a  choice  by  entering  the  associated  number.   This 
step  allows  the  user  to  choose  the  method  used  by  the 
program  to  obtain  discrete  frequency  values  used  to 
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generate  data  points.   LOGARITHMIC  INTERPOLATION  pro- 
duces data  points  that  appear  equally  spaced  on  a 
logarithmic  scale  between  the  minimum  and  maximum 
frequency  values  entered  in  steps  6  and  7.   With  this 
choice  a  Bode  and/or  Nyquist  plot  may  be  obtained. 
LINEAR  INTERPOLATION  produces  data  points  that  appear 
equally  spaces  on  a  linear  scale.   Only  a  Nyquist 
plot  may  be  obtained  with  this  choice.   DISCRETE  VALUES 
SUPPLIED  allows  the  user  to  supply  a  list  of 
discrete  frequency  values  of  interest.   No  plots 
may  be  obtained,  however,  with  this  option.   Below 
is  a  table  summarizing  the  available  graphical  output 
for  each  choice. 

CHOICE     BODE     NYQUIST 


0 

YES 

YES 

1 

NO 

NO 

2 

NO 

YES 

STEP  10  (assumes  choice  1  in  step  9): 

SCREEN  PROMPT: 

TYPE  IN  THE  DISCRETE  FREQUENCY  VALUES, 
(number  entered  in  step  8)  VALUES  NEEDED 
FREQ  (1)? 

REMARKS : 

Enter  the  list  of  frequency  values  as  prompted.   When 

the  last  value  has  been  entered  the  program  will  run 

to  completion  outputting  data  similar  to  Figure  4  with 

the  exception  of  the  plots. 
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STEP  11  (assumes  choice  0  in  step  9): 
SCREEN  PROMPT: 

BODE  PLOT?  (Y/N)  -> 
REMARKS : 

The  program  outputs  tabular  data  regardless  of  the 
choice  made  here. 
STEP  12  (assumes  choice  0  or  2  in  step  9): 
SCREEN  PROMPT: 

NYQUIST  PLOT?  (Y/N)  -> 
REMARKS : 

The  program  outputs  tabular  data  regardless  of  the 
choice  made  here. 

After  this  step,  the  program  runs  to  completion. 
See  sample  output  in  Figure  4. 
3.   FRESP  Example 

It  is  desired  to  know  the  frequency  response  of  the 
transfer  function 

8(0.5  +  s) 


G(s)  = 


4  +  6s  +  3s2  +  s3 


The  frequency  range  of  interest  is  GJ  =  0.1  rad/sec  to  100 
rad/sec.   A  Bode  plot  of  the  magnitude  and  phase  is  desired 
but  a  Nyquist  plot  is  not.   The  plots  should  be  generated 
from  100  values  of  frequency  logarithmically  spaced  from 
0.1  to  100. 


59 


Referring  to  the  steps  described  in  the  previous 
section,  FRESP  User's  Guide,  this  problem  would  be  entered 
as  follows: 

STEP  1.   Enter  "EXAMPLE   CR  "  where  <CR>  =  carriage 
return. 

STEP  2.   Enter  "8  <CR>." 

STEP  3.   Enter  "F,l  <CR>." 

STEP  4B.   Enter  ".5  <CR>  0  <CR>." 

STEP  5.   Enter  "P,3  <CR>  4  <CR>  6  <CR>  3  <CR>." 

STEP  6.   Enter  ".1  <CR>." 

STEP  7.   Enter  "100  <CR>." 

STEP  8.   Enter  "100  <CR>." 

STEP  9.   Enter  "0  <CR>." 

STEP  11.   Enter  "Y  <CR>." 

STEP  12.   Enter  "N  <CR> . " 

The  program  will  then  continue  to  completion  pro- 
ducing the  output  seen  in  Figure  4. 

The  output  begins  with  the  program  name  and  problem 
identification.   The  transfer  function  gain  is  given  followed 
by  the  transfer  function  numerator  and  denominator  each 
listed  as  coefficients  in  ascending  powers  of  s  and  then  as 
the  real  and  imaginary  parts  of  the  roots. 

The  second  page  of  output  is  tabular  data.   It 
consists  of  the  radian  frequency  and  the  transfer  function's 
corresponding  real  and  imaginary  parts,  magnitude,  and  phase 
in  radians  and  in  degrees.   Although  100  frequency  values 
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were  generated  and  their  associated  data  points  plotted, 
tabular  data  appears  only  for  every  other  data  point.   To 
reduce  unnecessary  output  and  for  formatting  purposes  the 
following  scheme  is  used  for  tabular  data. 


NUMBER  OF  FREQUENCIES  TABULAR  DATA 

REQUIRED  FROM  STEP  8  PRINTED  FOR 

0  to  50  every  freq. 

51  to  100  every  other  freq 

101  to  150  every  3rd  freq. 

151  to  200  every  4th  freq. 


The  next  two  pages  of  output  are  the  Bode  plots  for 
amplitude  and  phase.   The  plots  are  headed  with  the  problem 
identification  and  type  of  plot  and  below  each  plot  is 
sufficient  information  to  allow  proper  interpretation  of  the 
plotted  data.   Note  that  the  phase  angles  are  normalized  to 
always  remain  between  -180  and  +180  degrees. 

This  is  the  end  of  the  output  generated  from  the 
example  as  it  was  input  since  a  Nyquist  plot  was  not  desired 
A  Nyquist  plot  is  included  as  the  last  page  of  Figure  4  for 
the  sake  of  completeness. 
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FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  EXAMPLE 

GAIN=3 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

.5 

1 


NUMERATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

-.5  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

4 
6 
3 

1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

-1  -1.73285081 

-1  1.7329  50  3  1 

-1  Q 


Figure  4   FRESP  Output 
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Figure    4    (Cont.)   PROBLEM    IDENTIFICATION    -    EXAMPLE 
RADIAN  FREQ.       REAL  PART  IMAGINARY  PART     MAGNITUDE  PWSE(RAD)  PHASE(CEG) 


.187226722 

.123284474 

.141747414 

.142975884 

.187381743 

.215443449 

.247787434 

.284893588 

.327454917 

.374493582 

.43287413 

.497702359 

.572234749 

.457933228 

.754443332 

.849749008 

1.80008001 

1.149757 

1.32194114 

1.51991109 

1.74752841 

2.00923302 

2.31012972 

2.45408731 

3.85385554 

3.51119177 

4.0370173 

4.44158888 

5.33449929 

4.13590734 

7.85488239 

3.1113884 

9.32483358 

18.7224724 

12.3284474 

14.1747413 

14.2975884 

18.7381745 

21.5443472 

24.7707439 

28.4803591 

32.7454921 

37.4493587 

43.2874135 

49.7702345 

57.2234774 

45.7933234 

75.4443342 

84.9749019 

180.300002 


1.81785827 

1.82247187 

1.82954848 

1.83885884 

1.35894983 

1.34445859 

1.08498827 

1.11285044 

1.14575147 

1.18495955 

1.23774441 

1.29903747 

1.37188815 

1.45241093 

1.53955702 

1.42457894 

1.49230749 

1.71474725 

1.44344872 

1.4103491 

.945794058 

.330589559 

-.134534879 

-.423812809 

-.50737737 

-.473084414 

-.484413557 

-.32894124 

-.259188135 

-.281847424 

-.154424892 

-.118278988 

-.898141942 

-.3485723775 

-.3528722615 

-.8395884317 

-.829941485 

-.8224335834 

-.8171782383 

-.9138054894 

-9.84425437E-03 

-7. 458249 18E-03 

-5.43781714E-83 

-4.2459 1987E-03 

-3.22744824E-83 

-2.44194898E-03 

-1. 34746 144E-83 

-1.39745541E-03 

-1.85734283E-03 

-7.99879835E-04 


.8514283748 

.8584342913 

.8463192482 

.874488284 

.3832923195 

.9928537744 

.188388178 

.187134157 

.111131384 

.110214475 

.191419548 

.3384138318 

.8422474713 

-.9288852383 

-.  118848545 

-.248481144 

-.441538446 

-.732787826 

-1.87484971 

-1.45814418 

-1.7451883 

-1.7944947 

-1.5559544 

-1.14877489 

-.383954387 

-.52954775 

-.342379248 

-.229992452 

-.142544144 

-.9921935898 

-.85981842 

-.3339127642 

-.8253780175 

-.3145334577 

-.3108533543 

-7.1 1247 154E-03 

-4.44519419E-83 

-3.8422238E-83 

-2.81114273E-93 

- 1.32 1427 14E-93 

-8.48521474E-04 

-5.79984204E-04 

-3.75450085E-94 

-2.4491 1144E-84 

-1.42394323E-04 

-1.04818434E-04 

-7.02459939E-95 

-4.42235779E-05 

-3.84084457E-05 

-2.08851985E-95 


1.91834782 

1.82415889 

1.83178223 

1.94152594 

1.85424531 

1.37841542 

1.89152785 

1.11799565 

1.1511236 

1.19236571 

1.24189277 

1.39153689 

1.37165954 

1.45256198 

1.54497711 

1.6453479 

1.75411404 

1.36478998 

1.96392254 

2.92239964 

1.99459581 

1.83654746 

1.54174192 

1.24324255 

.'959479441 

.713394931 

.531384177 

.394282748 

.295723333 

.221178139 

.145791393 

.124515475 

.8934454855 

.9785491932 

.9531914145 

.349135442 

.8383024494 

.8223892471 

.8172954174 

.0130724492 

9. 8824953 1E-03 

7.4721 1725E-83 

5.45838488E-83 

4.2738587E-83 

3.23173108E-03 

2. 444294 15E-93 

1.84879722E-03 

1.39841954E-03 

1. 9577380 1E-93 

8.80129944E-84 


.0507189222 

.8572328479 

.8443257541 

.3714873348 

.3798898087 

.0848884134 

.8920249474 

.3959743558 

.896691838 

.8925996297 

.3817563658 

.9619339349 

.9383193662 

-.8143737124 

-.8765349743 

-.159956477 

-.264252852 

-.483358573 

-.579143334 

-.799389724 

-1.84534155 

-1.34239719 

-1.45784784 

-1.91844312 

-2.13377179 

-2.38584541 

-2.44182781 

-2.55882432 

-2.43842774 

-2.71143854 

-2.77251189 

-2.32375481 

-2.34722885 

-2.98430799 

-2.9340947 

-2.94344184 

-2.98782444 

-3.8874041 

-3.3250474 

-3.04033484 

-3.35359432 

-3.06510255 

-3.07539595 

-3.88377726 

-3.89132896 

-3.39737757 

-3.18357722 

-3.18853252 

-3.11284111 

-3.11658766 


2.99593123 

3.2320646 

3.68559544 

4.1373332 

4.53146882 

4.93258435 

5.27275758 

5.4993249 

5.54883621 

5.3859342 

4.48429633 

3.55124725 

1.74584839 

-.323339329 

-4.33513258 

-9.11326839 

-15. 2551243 

-23,1389419 

-33.1325894 

-45.7979991 

-61.3395965 

-78.3424434 

-94.9418367 

-139.931333 

-122,256162 

-132.379573 

-139.363595 

-146.135483 

-151.132287 

-155.3655 

-158.353237 

-161.739291 

-144.279447 

-144.484593 

-148.226809 

-169.792724 

-171.143947 

-172.311739 

-173.322511 

-174.193418 

-174.958129 

-175.417593 

-174.198933 

-174.487435 

-177.119788 

-177.495374 

-177.321939 

-173.352722 
-173.547383 


Figure   4 
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Figure  4  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  BODE  PLOT  (AMPLITUDE)  XX 


l    l   l   l  l  J  l  l  l 


T 1  I  I  I  I  I  I 


T I 1  I  I  I  I 


i 


'   '  '  I  '  '  U 


_L 


J i '  '  '  I  I  ' 


J I I  I  I 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  = 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  = 


.  1  RADIANS/SEC 
198  RADIANS/SEC 


AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-80  DECIBELS 


Figure  4 
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Figure  4  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  BODE  PLOT  (PHASE)  XX 


T 1 1  MINI 


T 1 1  I  1  I  I  I  I 


T 1 1  I  I  I  I  I 


JL 


i 


J '  '  LLP 


J I    I  I  I  I  I 


T1— +■ 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  90 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  =  90 


DEGREES 

=  . 1  RADIANS/SEC 
=  180  RADIANS/SEC 
DEGREES 


MINIMUM  PHASE  ON  ORDINATE  SCALE  =  -180  DEGREES 


Figure  4 
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Figure  4  (Cont . ) 


PROBLEM  I DENTIFI CATION  -  EXAMPLE 

XX  NYQUIST  PLOT  XX 


ABSCISSA  ->  REAL  PART  OF  6<JN) 
ORDINATE  ->  IMAGINARY  PART  OF  G<JW) 
TIC  MARKS  SHOW  INTERVALS  OF  UNITY 
AXES  CROSS  AT  ORIGIN 


Figure  4 
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B .   RTLOC 

1.   Root  Locus  Program  (RTLOC) 

a.  Introduction 

When  analyzing  a  system  it  is  helpful  to  know  the 
location  of  the  closed  loop  poles  in  the  s-plane.   These 
closed  loop  poles  are  the  roots  of  the  characteristic 
equation  and  determine  the  basic  characteristics  of  the 
transient  response  of  a  closed  loop  system.   A  root  locus 
plot  is  a  plot  of  the  roots  of  the  characteristic  equation, 
usually  as  a  function  of  the  gain  of  the  transfer  function, 
and  therefore  it  is  a  valuable  tool  for  system  analysis. 

b.  Description  of  Program 

The  RTLOC  program  [Ref.  1:   pp.  114-121]  calcu- 
lates and  plots  the  roots  of  the  equation 

1  +  KG(s)  =  0 
as  a  function  of  K.   G(s)  is  assumed  to  be  a  rational  function 
of  the  form 

G(s)  =   N(S) 


D(s) 


and  the  root  locus  becomes  the  locus  of  roots  of  D„(s)  as  K 

1\ 


varies  where 


D  (s)  =  D(s)  +  KN(s) 


RTLOC  uses  an  algebraic  plus  linear  progression 
scheme  to  vary  K  to  give  reasonable  spacing  of  the  roots. 
D„(s)  is  obtained  for  each  value  of  K  and  the  subroutine  PROOT 
is  used  to  calculate  its  roots. 
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The  scheme  used  to  calculate  values  of  K  assumes 
that  K  takes  on  only  positive  values.   If  K  is  to  range 
through  negative  value,  the  value  of  smaller  magnitude  (less 
negative)  must  be  used  as  the  minimum  value.   The  routine 
starts  at  the  maximum  (less  negative)  value  and  becomes 
increasingly  negative  until  the  lower  limit  is  reached. 
If  both  positive  and  negative  values  are  desired  then  two 
separate  runs  must  be  made  with  only  positive  and  only 
negative  values. 

The  RTLOC  program  has  a  feature  that  allows  the 
user  to  specify  a  range  of  s  and  (jj    values  around  an  area 
of  the  root  locus  plot  that  is  of  interest.   That  area  is 
then  enlarged  and  more  closely  spaced  values  of  K  are  generated 
giving  more  detail  to  the  plot. 

The  major  subroutines  used  are  PROOT  and  SEMBL. 
The  subroutine  SPLIT  used  in  the  FORTRAN  version  is  replaced 
by  additional  coding  in  the  main  program  and  by  various 
plotting  subroutines. 

PROOT.   This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.   It  is  used  to  determine 
the  roots  of  the  numerator  and  denominator  of  G(s)  when  it 
is  entered  in  polynomial  form.   It  is  also  used  to  determine 
the  roots  of  D„(s)  for  each  value  of  K  generated. 

SEMBL.   This  subroutine  determines  the  coeffi- 
cients of  a  polynomial  from  the  roots  of  the  polynomial.   It 
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is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  G(s)  when  they  are  entered  in 
the  factored  form.   This  subroutine  together  with  PROOT 
provides  the  feature  that  G(s)  may  be  entered  in  either 
factored  or  coefficient  form.   It  appears  in  the  output  in 
both  forms. 

c.   Program  Translation  Problems 

In  addition  to  the  programing  considerations 
discussed  in  section  II,  two  problems  were  encountered. 
When  a  range  of  negative  gains  is  entered  the  program  cal- 
culates only  one  value  of  K  due  to  the  logic  of  line  229  in 
the  FORTRAN  code.   This  works  correctly  only  for  positive 
ranges  of  gain.   (Note  that  K  is  represented  by  the  variable 
G  in  the  computer  program. )   In  the  BASIC  version  coding 
was  added  to  test  the  sign  of  the  range  of  gains  and  modify 
the  logic  to  correctly  handle  the  negative  case. 

The  second  problem  was  encountered  when  using 
gains  of  very  small  magnitude  necessary  in  using  this 
program  for  w1 -plane  analysis.   The  schemes  used  to  generate 
values  of  K  for  a  reasonable  spacing  of  the  roots  when  plotted 
worked  fine  for  the  magnitudes  usually  encountered  in  s-plane 
analysis  but  was  not  flexible  enough  to  handle  gains  of  much 
smaller  magnitudes  (e.g.  0  to  1E-3).   To  add  this  necessary 
flexibility  the  following  FORTRAN  lines  used  to  generate 
values  of  gain  G, 


69 


227  G=1.15*(G+SIGNG*0.05) 

228  G=1.04*(G+SIGNG*0.02) 

were  modified  by  replacing  the  constant  values  0.05  and  0.02 

with  the  variables  Dl  and  D2 ,  respectively,  where 

D1=ABS ( GMIN-GMAX ) / 700 
D2=ABS(GMIN-GMAX)/1500 

so  that  the  gain  increment  was  a  function  of  the  range  of 

gains  of  interest . 

2.   RTLOC  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

ROOT  LOCUS  PROGRAM  (RTLOC) 

THIS  PROGRAM  PLOTS  THE  ROOT  LOCUS  OF  A 
DESCRIBED  SYSTEM. 

PROBLEM  IDENTIFICATION  -> 
REMARKS : 

Label  the  problem.   Type  in  any  appropriate  combina- 
tion of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.   This  input  is  limited  to  255 
characters . 
STEP  2: 

SCREEN  PROMPT: 

INPUT  THE  RANGE  OF  GAINS  (MIN,MAX)  -> 
REMARKS : 

Enter  the  range  of  transfer  function  gains  of  in- 
terest.  Note  that  the  minimum  gain  is  the  gain  of 

lowest  absolute  value. 
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STEP  3: 

SCREEN  PROMPT: 

•KEY1  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
'KEY1  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY1,  ORDER  FOR  NUMERATOR,  N(S)  -> 

REMARKS : 

Choose  the  preferred  method  of  entering  the 

numerator  polynomial  of  the  transfer  function 

G(s). 

STEP  4A  (assumes  coefficient  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  COEFFICIENTS 

INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  s  OR 

w1 

INPUT  COEFFICIENT  OF  S~0  -> 

REMARKS : 

Enter  the  coefficient  along  with  its  algebraic  sign 

as  prompted.   The  program  assumes  that  the 

coefficient  of  the  highest  power  of  s  is  one. 

STEP  4B  (assumes  factored  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 

FOR  FACTOR  1 
REAL  PART  -> 

IMAGINARY  PART  -> 
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REMARKS : 

Enter  the  factor  with  its  algebraic  sign.   Do  not 

enter  the  associated  root.   (e.g.  For  "(s-1)"  enter 

"-1".)   After  the  real  part  of  the  factor  is 

entered  the  program  asks  for  the  imaginary  part 

allowing  for  input  of  quadradic  factors.   If  a 

non-zero  imaginary  part  is  entered  the  program 

automatically  enters  its  complex  conjugate. 

STEP  5: 

SCREEN  PROMPT: 

•KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR,  D(S)  "> 

REMARKS : 

The  denominator  is  entered  just  as  described  for  the 

numerator  in  steps  3  and  4. 

STEP  6: 

SCREEN  PROMPT: 

WOULD  YOU  LIKE  TO  LOOK  AT  ONLY  A  PART  OF  THE 
ROOT  LOCUS?   (Y/N)  -> 

REMARKS : 

Enter  an  "N"  for  the  option  of  viewing  the  root  locus 

for  the  entire  range  of  gains  entered  in  STEP  2. 

Enter  a  "Y"  for  the  option  of  viewing  only  a  portion 

of  the  root  locus  defined  by  minimum  and  maximum 

values  of  sigma  and  (jj   of  interest  to  the  user. 

Usually,  it  is  best  to  view  the  entire  root  locus 
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first  then  decide  if  any  section  needs  to  be 
enlarged  to  reveal  more  detail  and  rerun  RTLOC 
to  view  only  that  portion. 
STEP  7  (assumes  "YM  entered  in  step  6): 
SCREEN  PROMPT: 

ENTER  SIGMA  MIN,  SIGMA  MAX  -> 
REMARKS : 

Enter  the  minimum  and  maximum  values  of  sigma  (i.e. 
the  real  axis)  of  the  portion  of  the  root  locus 
plot  of  interest. 
STEP  8  (assumes  "Y"  entered  in  step  6): 
SCREEN  PROMPT: 

ENTER  OMEGA  MIN,  OMEGA  MAX  -> 
REMARKS : 

Enter  the  minimum  and  maximum  values  of  omega  (i.e. 
the  imaginary  axis)  of  the  portion  of  the  root 
locus  plot  of  interest. 
STEP  9: 

SCREEN  PROMPT: 

PRINT  OUT  OF  GAIN  DATA?  (Y/N)  -> 
REMARKS : 

This  gives  the  user  the  option  to  suppress  the  print 
out  of  the  gain  data  if  it  is  not  needed.   This  saves 
time  and  paper  if  only  the  root  locus  plot  is  desired 
After  this  step  the  program  runs  to  completion. 
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3.   RTLOC  Example 

It  is  desired  to  know  the  root  locus  for  the  open 
loop  transfer  function  G(s)  given  by 


G(s)  =  1^_jL_s 


s(8  +  9s  +  s2) 


for  gain  K  of  from  0  to  30.   Referring  to  the  steps  des- 
cribed in  section  2,  RTLOC  User's  Guide,  this  problem  would 
be  entered  as  follows: 

STEP  1.   Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 
return . 

STEP  2.   Enter  "0,30  <CR>." 

STEP  3.   Enter  "F,l  <CR>." 

STEP  4B.   Enter  "1.2  <CR>  0  <CR> . " 

STEP  5.   Enter  "F,3  <CR>  0  <CR>  0  <CR>  1  <CR>  0  <CR> 
8  <CR>  0  <CR>." 

STEP  6.   Enter  "N  <CR>." 

STEP  9.   Enter  "Y  <CR>." 

The  program  will  then  run  to  completion  producing 
the  output  seen  in  Figure  5. 

The  output  begins  with  the  program  name  and  problem 

identification.   Next  the  transfer  function  numerator  and 

denominator  are  listed  as  coefficients  in  ascending  powers 

of  s  and  then  as  the  real  and  imaginary  parts  of  the  roots 

(i.e.  open  loop  zeros  for  the  numerator  and  open  loop  poles 

for  the  denominator).   The  minimum  and  maximum  values  of  gain 

are  listed  next. 
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If  the  option  to  view  only  a  part  of  the  root  locus 
is  selected,  this  fact  is  stated  next  followed  by  the  ranges 
of  sigma  and  omega  that  are  chosen.   (See  Figure  6) 

The  next  several  pages  of  output  contain  the  gain 
data  which  may  be  suppressed  with  the  appropriate  response  in 
step  9.   The  data  point  number  is  listed  with  the  value  of 
gain  and  the  real  and  imaginary  parts  of  the  corresponding 
roots  of  the  open  loop  system. 

The  last  page  of  output  contains  the  root  locus  plot 
itself.   It  begins  with  the  problem  identification  and 
heading  and  ends  with  sufficient  data  listed  to  allow  proper 
interpretation  of  the  plot. 

Figure  6  shows  the  output  from  the  same  example 
problem  with  the  gain  data  suppressed  and  the  option  to 
view  only  a  part  of  the  root  locus  chosen. 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EXAMPLE 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

1.2 
1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-1.2  0 


DENOMINATOR    COEFFICIENTS    -    IN    ASCENDING    POWERS    OF    S    OR    W' 

9 
3 
9 

1 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 

-1  9 

-S  0 

0  0 

MIN  GAIN  MAX  GAIN 

0  30 


Figure  5   RTLOC  Output 
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Figure  5  (Cont . ) 


1 


GAIN  =  9 


3 


18 


GAIN  =  .9  492857  143 


GAIN  =  . 19  5964236 


GAIN  =  .17  1  144643 


GAIN  =  .24619  2954 


GAIN  =  .33239  39  76 


GAIN  =  .431434251 


GAIN  =  .545435193 


GAIN  =  .676536933 


GAIN  =  .32739221 


ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1  9 

-3  9 

9  9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-7.49  391563E-9  3  9 
-.998530111  8 
-7.9940  1093        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.0159695125  9 
-.99691764  1  8 
-7.93711234        8 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.9  258639696  9 
-.994965313  8 
-7.97916572        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.937328949  0 
-.992664433  9 
-7.97999752        9 

ROOTS  ARE 
REAL  PART  I MAG.  Ph RT 

-.0506034588  0 
-.939941249  8 
-7.95945829        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.8668224557  3 
-.98678  1776  0 
-7.94727577        8 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.0339456946  9 
-.9S2324229  0 
-7.9332300  8        8 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.104334979  0 
-.973147608  0 
-7.91701742        8 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.972453345  8 
-7.89S29233  8 
-.129253772        8 


Figure  5 
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Figure  5  (Cont . ) 
1 1     GAIN  =  1 .90868326 


12     GAIN  =  1 .20907146 


13     GAIN  =  1 .4293678? 


14     GAIN  =  1.6930587? 


15     GAIN  =  1  .9963033: 


16 


17 


13 


19 


GAIN  =  2.34503453 


GAIN  =  2.7460  7543 


GAIN  =  3.20727245 


GAIN  =  3.73764904 


20     GAIN  =  4.34753211 


ROOTS  ARE 


REAL  PART 

I  MAG. 

PART 

-.965438583 

8 

-7.37665069 

0 

-.  15791073 

8 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.956660582 

0 

-7.3516177 

0 

-.  191721725 

8 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.945447526 

8 

-7.82263429 

0 

-.231913189 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.7390  39  12 

8 

-.230258299 

0 

-.930702531 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.75004614 

0 

-.339507135 

8 

-.910446682 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.70471555 

0 

-.414324269 

8 

-.330460  137 

8 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.65191587 

0 

-.520147285 

0 

-.327936925 

8 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.59027312 

0 

-.704863441 

-. 10 1 133243 

-.704863441 

.10  11 

33243 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.51810122 

0 

-.740949392 

-.213 

123799 

-.740949392 

.2131 

23799 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-7.43330  403 

0 

-.783347989 

-.297 

0  19334 

-.733347989 

.2978  19334 

Fi 

gure  5 
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Figure  5  (Cont .  ) 

21     GAIN  =  5.84908514 


22     GAIN  =  5.35564162 


!3     GAIN  =  6.7332735S 


24     GAIN  =  7.35995833 


25     GAIN  =  9.0  763435? 


26     GAIN  =  10.4876559 


27     GAIN  =  12. 1 190899 


!8     GAIN  =  13.975839  1 


29     GAIN  =  16. 1215532 


30     GAIN  =  18.5390777 


ROOTS  ARE 
REAL  PART  I MAG.  PART 

-7.33323262  9 
-.333333692        -.362332638 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-7.2144639         9 
-.892765554        -.420657492 
-.392765554        .420657402 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-7.07243055        9 
-.963759725        -.471271626 
-.963759725        .47  1271626 

ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-6.90  10  354         9 
-1.0494823         -.51342799  3 
-1.0494323         .513427903 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-6.69112225        0 
-1.15443333        -.54325973  1 
-1.15443333        .543259731 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-6.42372313        9 
-1.23563594        -.5529  76334 
-1.23563594        .552976334 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-6.9  3944969        9 
-1.45527515        -.513289319 
-1.45527515        .513289  319 

ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-5.62140  192        9 
-1.68929994        -.369  139954 
-1.63929994        .369  139054 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-2.6433642  9 
-4.34656548  9 
-1.51007033        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-3.30  14976  1        -1.23137034 
-3.80  14976  1        1.231370  34 
-1.39700479        9 

Figure  5 
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Figure  5  (Cont . ) 

31     GAIN  =  21.426725 


32     GAIN  =  24.6990  195 


33     GAIN  =  23.4428031 


ROOTS  ARE 
REAL  PART 
-3.82838139 
-3.82338189 
-1.34323622        9 

ROOTS  ARE 


I MAG.  PART 


.  1  1737932 


REAL  PART 
-3.8449632 
-3.3449632 
-1 .3100636 


I MAG.  PART 
-2.79355914 
2.79855914 
9 


ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-3.3563616  1        -3.4  1216555 
-3.3563616  1        3.4  1216555 
-1.23727673        0 


Figure  5 
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Figure  5  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  ROOT  LOCUS  PLOT  XX 


1        1 

i        T        r 

r       r 

m 

♦ 

♦ 

- 

♦ 

- 

♦ 

♦~s 

- 

*U 

♦ 

- 

♦ 

- 

- 

- 

.1 1 

!                 1                 I 

■  i i 

ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -8  TO  1 

ORDINATE,  -4  TO  5 


Figure  5 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EXAMPLE  WITH  OPTION 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

1.2 

1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-1.2  9 


DENOMINATOR    COEFFICIENTS    -    IN    ASCENDING    POWERS    OF    S    OR    W 

9 
S 
? 

1 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 


-1  9 

-8  9 

9  9 


MIN  GAIN  MAX  GAIN 

9  39 


^X^X^^XXaxXX^XXX^^XM^^^^XXKX*^^-^-(x^^^-(¥¥x^^^x^***S 


OPTION  TAKEN 

SIGMA  MIN  =  -2  SIGMA  MAX  =  9 

OMEGA  MIN  =  -1  OMEGA  MAX  =  1 

Figure  6   RTLOC  Output  with  Option 
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Figure  6  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EXAMPLE  HITH  OPTION 
XX  ROOT  LOCUS  PLOT  XX 


-L 


ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -2  TO  1 

ORDINATE,  -1  TO  2 


Figure  6 
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V.   w1 -PLANE  ANALYSIS 

Digital  control  systems  are  becoming  more  and  more  common. 
Digital  control  laws  have  unique  characteristics  that  can 
only  be  approximated  by  using  classical  techniques  in  the 
continuous  s  domain. 

In  the  w*  domain  all  analog  control  system  design  tech- 
nology transfers  completely  for  digital  control  system  design. 
An  important  advantage  of  the  w1  domain  is  that  non-minimum 
phase  effects  of  the  sampling  and  data-hold  operations  and 
of  sampling  rate  can  be  directly  accounted  for  without 
approximation  while  using  conventional  frequency  domain 
design  and  analysis  tools  such  as  root  locus  and  Bode 
plots.   [Ref.  4] 
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A .   BACKGROUND 

This  s  domain  is  used  for  continuous  system  analysis. 
When  a  digital  system  is  considered  the  z,  w,  or  w1 domain 
must  be  used.   These  domains  are  related  as  shown  below. 

sT 
z  =  e 


z-1        sT 

w  =  =  tanh  — 

z+1        2 


w'  =  (2/T)  w 
where  T  is  the  sampling  period  and 


tanh  x  = 


x    x 
e   -  e 


x  ,   x 
e   +  e 


In  an  s-plane  root  locus  plot  the  region  of  stability  is 
the  left  half  plane,  that  is  all  roots  with  negative  real 
parts.   This  region  of  stability  is  mapped  into  a  unit  circle 
with  its  center  at  the  origin  in  the  z-plane.   3y  use  of 
the  bilateral  transformation  shown  above,  the  stability  region 
of  the  z-plane  is  mapped  back  into  the  left  half  plane  to 
form  the  w-plane.   The  w! -plane  takes  it  a  step  farther  by 
multiplying  the  w-plane  by  the  factor  2/T  where  T  is  the 
sampling  period.   This  gives  the  w1 -plane  the  property 
that  not  only  is  the  region  of  stability  the  left  half  plane 
as  in  the  continuous  s-plane  but  w1  approaches  s  as  the  sam- 
pling period  T  approaches  zero  [Ref.  5]. 

85 


B.   TRANSFER  FUNCTION  APPROACH  USING  THE  w1 -PLANE 
In  this  section  the  second  order  system 


>  K 


-> 


10 


s(s+10) 


is  analyzed  using  the  two  transfer  function  programs  dis- 
cussed in  part  IV  of  this  thesis.   It  is  then  converted 
to  the  w'-plane  and  analyzed  for  periods  of  .001  seconds, 
.01  seconds,  and  .1  seconds  using  the  same  two  transfer  func- 
tion programs.   These  results  are  used  to  gain  insight  into 
w'-plane  analysis. 

The  system  is  converted  to  the  w'-plane  by  adding  a 
sampler  and  a  digital  to  analog  converter  in  the  form  of  a 
zero  order  hold.   This  modified  system  is  shown  below. 


7\       **         v 

-sT 
1-e 

s 

s 

10K 

v. 

>o— — r             •? 

s 

s(s+10) 

The  new  open  loop  transfer  function  is  transformed  to 
the  z-plane  and  then  to  the  w'-plane  as  follows: 


G(s)  =  10K(l-e~sT) 


s  (s+10) 
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G(z)  =  10K(l-z  1) 


10 


Tz 


,-   -10T\ 
(1-e     )z 


(z-1)2    10(z-l)(z-e  10T) 


=  K 


1-e 


-10T 


.,    -rt/    -10T\ 
z-1    10(z-e     ) 


For  the  period 


T=.001 


G(z)  reduces  to 


G(z)  = 


4.98337  x  10  6  K(z+. 996681) 
(z-l)(z-. 990049834) 


(1) 


Similarly,  for 


T=.01 


G(z)  = 


4.837418  x  10~4  K( z+. 9672185 ) 


(z-l)(z-. 904837418) 


And  for 


T=.l 


G(z)  = 


.03678794412  K( z+. 718281827) 
(z-l)(z-. 3678794412) 


Now  G(z)  is  converted  to  the  w'  domain  using  the  rela- 
tionship 
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1  +  T/2  W 

z  =  

1  -  T/2  w1 

For  period 

T=.001 


1  +  .0005  w1 
z  =  


1  -  .0005  w* 


Substituting  into  equation  (1)  gives 


/  l+.0005w' 

4.98337  x  10    K    ■ 

V  l-.0005w' 
G(w')  = 

l+.0005w'     \  /l+.0005w' 

-  1     • 

1-.0005W'     /  V  l-.0005w' 


which  reduces  to 


-4.156878  x  10    K  (wf  +  1202819)  (w1  -  2000) 

G(w')  =  (2) 

w'  (w«  +  9.9999) 


Note  that  the  gain  is  1.0K  which  is  the  same  as  the  s  domain 
gain.   Similarly,  for 

T=.01 


l+.005w' 

z=  


1-.005W' 
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-4.1625027  x  10    K(w' +12002 .00418 ) (w ' -200 ) 

G(w')  =  (3) 

w»  (w1  +  9.991674985) 


And  for 


T=.l 


l+.05w» 
z  =  


l-.05w 


-.0037882843  K(w'  +  121 . 9858708 ) (w' -20 ) 
G(w')  =  (4) 

w'  (w'  +  9.242343139) 


Note  again  that  the  gain  is  1.0K  also  in  equations  (3)  and 
(4). 

1.   Frequency  Response 
a.   s-Plane 

The  open  loop  transfer  function  of  the  system 
in  the  s  domain  is 


10K 
G(s)  =  


s(s+10) 


The  frequency  response  of  G(s)  with  a  gain  K  of  5  is  described 
by  the  output  of  the  FRESP  program  shown  in  Figure  7.   From 
the  Bode  plots  it  is  seen  that  the  gain  margin  is  infinite 
and  the  phase  margin  is  65.1  degrees. 
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b.  w'-Plane  with  a  Period  of  .001  Seconds 

The  open  loop  transfer  function  for  the  equiva- 
lent sampled  system  is  represented  by  equation  (2)  for  a 
period  T  of  .001  seconds.   This  transfer  function  is 
entered  into  the  FRESP  program  as  if  w1  were  an  s.   The 
resulting  output  in  Figure  8  is  interpreted  as  described 
earlier  in  the  thesis  keeping  in  mind  that  all  references 
to  s  are  actually  referring  to  w1.   From  this  output  it 
is  seen  that  the  gain  margin  is  no  longer  infinite  as  it 
was  in  the  continuous  case  but  it  is  still  quite  high  at  50 
dB.   The  phase  margin  has  dropped  slightly  from  65.1  degrees 
to  64.6  degrees. 

c.  w'-Plane  with  a  Period  of  .01  Seconds 

If  the  sampling  rate  of  the  system  is  decreased 
so  that  the  sampling  period  is  .01  seconds,  then  the  open 
loop  transfer  function  is  represented  by  equation  (3). 
Inputting  the  transfer  function  into  FRESP  results  in  the 
output  shown  in  Figure  9.   From  the  Bode  plots  the  gain 
margin  is  found  to  be  32  dB.   This  is  down  from  50  dB  for 
the  case  of  a  period  of  .001  seconds  indicating  a  decrease 
in  stability.   The  phase  margin  is  64.5  degrees  which  is  also 
lower  than  the  previous  case  but  only  by  one  tenth  of  a 
degree. 

d.  w'-Plane  with  a  Period  of  .1  Seconds 

If  the  period  is  again  increased  to  .1  seconds, 
the  open  loop  transfer  function  is  represented  by  equation 
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(4).   Entering  this  transfer  function  into  FRESP  produces 

the  output  seen  in  Figure  10.   From  this  output  the  gain 

margin  is  found  to  have  decreased  to  13.2  dB  and  the  phase 

margin  to  53.7  degrees. 

e.   Summary  of  Frequency  Response  Results 

The  table  below  is  a  brief  summary  of  the  gain 

and  phase  margins  found  in  each  case. 

CASE  GAIN  MARGIN        PHASE  MARGIN 

s-plane  infinite  65.1  deg. 

w1 ,  T=.001  50  dB  64.6  deg. 

w1 ,  T=.01  32  dB  64.5  deg. 

w1 ,  T=.l  13.2  dB  53.7  deg. 

It  can  be  seen  from  this  table  that  the  con- 
tinuous case  is  the  most  stable.   The  sampled  cases  become 
less  stable  as  the  sampling  period  increases.   It  is  also 
noticed  that  the  gain  margin  is  more  sensitive  to  changes 
in  the  sampling  period  than  the  phase  margin. 
2 .   Root  Locus 

Now  a  comparison  is  made  using  the  root  locus  program 
using  the  same  cases  used  above, 
a.   s-Plane 

If  the  open  loop  transfer  function 


10K 
G(s)  =  


s(s+10) 


is  entered  into  the  RTLOC  program,  the  result  is  the  output 
in  Figure  11.   The  root  locus  plot  shows  that  the  system 
never  becomes  unstable  at  any  gain. 
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b.  w* -Plane  with  a  Period  of  .001  Seconds 

If  the  equivalent  sampled  transfer  function  with  a 
sampling  period  of  .001  seconds  represented  by  equation  (2) 
is  entered  into  RTLOC ,  the  output  will  be  that  found  in 
Figure  12.   From  the  root  locus  plot  it  can  be  seen  that 
there  is  a  slight  tendency  for  the  plot  to  curve  toward  the 
right  half  plane  as  it  moves  further  from  the  real  axis.   The 
tendency  is  so  slight  that  it  will  take  a  very  large  gain 
to  drive  the  system  unstable. 

c.  w' -Plane  with  a  Period  of  .01  Seconds 

If  the  sampling  period  is  increased  to  .01  seconds 
the  transfer  function  is  represented  by  equation  (3). 
Entering  this  transfer  function  into  RTLOC  results  in  the 
output  seen  in  Figure  13.   It  can  be  seen  that  the  tendency 
for  the  curve  to  bend  toward  the  unstable  right  half  plane 
is  increased  indicating  that  a  lower  value  of  gain  than  in 
the  previous  case  will  drive  the  system  to  instability. 

d.  w'-Plane  with  a  Period  of  .1  Seconds 

When  the  sampling  period  is  further  increased  to  .1 
seconds  represented  by  equation  (4),  the  resulting  RTLOC 
output  is  that  seen  in  Figure  14.   In  the  root  locus  plot 
for  this  case  the  tendency  to  become  unstable  is  much  more 
pronounced.   Here  the  plot  bends  into  the  right  half  plane 
within  the  limited  boundaries  of  the  portion  plotted.   It 
becomes  unstable  for  values  of  gain  K  greater  then  24.6. 
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e.   Summary  of  Root  Locus  Comparison 

It  can  be  seen  that  the  continuous  system 

10K 
G(s)  =  


s(s+10) 


is  always  stable.   When  this  system  is  sampled  it  can  be 
seen  by  equations  (2),  (3),  and  (4)  that  a  zero  is  added 
in  the  right  half  plane  making  it  a  non-minimum  phase 
system.   It  can  also  be  seen  that  by  decreasing  or  increasing 
the  sampling  period  the  distance  of  this  zero  from  the 
origin  increases  or  decreases  respectively.   The  closer 
this  zero  is  to  the  imaginary  axis,  the  greater  effect  it 
has  on  bending  the  root  locus  into  the  right  half  plane. 
This  effect  can  be  seen  by  examining  the  root  locus  plots 
in  Figures  11  through  14. 
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FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  EX  1  S-PLANE 
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Figure  7   s-Plane  Frequency  Response  Example 
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Figure    7    (Cont . ) 

s  'PROBLEM    IDENTIFICATION    -    EX1    S-PLANE 
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Figure  7  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EX1  S-PLANE 
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Figure    7    (Cont .  ) 


PROBLEM    I DENTIFI CATION    -    EX  1    S-PLANE 
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Figure  7 
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Figure  8  w1 -Plane  Frequency  Response  Example,  T=.001 
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Figure  8  (Cont .  ) 
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Figure  8  (Cont .  ) 
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REAL  PART  IMAGINARY  PART 

-12002.0042  9 

200  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

0 

9.  991 67493 
1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

0  9 

-9.99167493  9 


Figure  9  w1 -Plane  Frequency  Response  Example,  T=.01 
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Figure    9    (Cont .  ) 


PROBLEM    I  DENT  I F I  CAT I  ON    -    EX  1    W ' -PLANE    T= . 9  1 


RADIAN  FREQ. 


REAL  PART 


IMGIMRY  PART     MGN1TUDE 


PHASE(RAD) 


PHASE (C£G) 


.187226722 

.123284474 

.141747414 

.162975084 

.187381743 

.21544346? 

.247797636 

.284893588 

.327454917 

.376493582 

.43287613 

.497782359 

.572236769 

.457933223 

.756463332 

.869749888 

1.88098801 

1.149757 

1.32194116 

1.51991109 

1.74752841 

2.88923302 

2.31812972 

2.65688781 

3.85385554 

3.51119177 

4.3378173 

4.64158888 

5.33449929 

6.13598734 

7.35488239 

8.1113084 

9.32603358 

10.7226724 

12.3284676 

14.1747418 

16.2975086 

18.7381745 

21.5443472 

24.7707639 

28.4883591 

32.7454921 

37.6493587 

43.2876135 

49.7702345 

57.2234774 

45.7933234 

75.6463342 

84.9749019 

100.900002 


-.524939547 

-.524920088 

-.524894344 

-.524860366 

-.524815428 

-.524756033 

-.524677539 

-.52457381 

-.524436749 

-.524255672 

-.524016491 

-.523780643 

-.523283494 

-.52273353 

-.522088018 

-.52105202 

-.519793484 

-.51813934 

-.515968635 

-.513126826 

-.509417817 

-.504596244 

-.49834874 

-.498350539 

-.488148485 

-.467296883 

-.451325398 

-.431817202 

-.40847643? 

-.381235751 

-.358349774 

-.316458228 

-.230573464 

-.244007236 

-.288143963 

-.174282984 

-.143437538 

-.116242624 

-.8929485537 

-.8734841787 

-.8575538451 

-.9447353616 

-.8345626152 

-.8265759352 

-.8283593365 

-.0155527033 

-.8118553889 

-9.82265597E-03 

-6. 8589325 1E-83 

-5.21006089E-03 


-46.4245341 

-40.5500649 

-35.2645454 

-30.6789756 

-26.6734542 

-23.1966296 

-28.1728793 

-17.5418867 

-15.2528989 

-13.2686853 

-11.527947 

-18.8288794 

-8.79747398 

-7.56513572 

-6.57018639 

-5.70343055 

-4.9479794 

-4.2891243 

-3.71405452 

-3.21161357 

-2.77289186 

-2.38784433 

-2.84914171 

-1.75212349 

-1.49852834 

-1.25981191 

-1.35419372 

-.876623211 

-.718746596 

-.588778561 

-.461388989 

-.359537358 

-.274264592 

-.28446584 

-.148767675 

-.185522173 

-.0728478818 

-.8488758445 

-.8317858744 

-.8197243394 

-.8115449147 

-6.14914418E-93 

-2.64835976E-83 

-4.59917689E-04 

8.48806795E-04 

1.57682656E-83 

1.9329344E-93 

2.85532664E-93 

2.9361368E-33 

1.93572135E-93 


46.6274892 

49.5534423 

35.2794713 

38.6754662 

26.6788168 

23.2825644 

20.1789016 

17.5488488 

15.2611045 

13.2718444 

11.5398588 

10.3337559 

8.72338387 

7.58317406 

6.59889981 

5.72718215 

4.97520709 

4.32930737 

3.74972327 

3.25234698 

2.31858943 

2.43979662 

2.1888924 

1.81944588 

1.56595572 

1.34368586 

1.14858184 

.97721191 

.826718149 

.694719472 

.579324936 

.47897969 

.392368887 

.313348855 

.255843176 

.283738772 

.  168884863 

.126899635 

.8932874138 

.0760854428 

.0587946732 

.9451564981 

.834663932 

.0265799146 

.9293749894 

.9156324331 

.9129119398 

9.2537932E-93 

7.15478944E-93 

5.55893487E-93 


-1.58205477 

-1.58374064 

-1.5854789 

-1.58790731 

-1.59944924 

-1.5934144 

-1.5968986 

-1.69969393 

-1.69516741 

-1.61931938 

-1.41622129 

-1.62391398 

-1.63981873 

-1.63978443 

-1.65998935 

-1.66199113 

-1.47544416 

-1.69191639 

-1.78883614 

-1.72922959 

-1.75253532 

-1.77911839 

-1.389367 

-1.34367592 

-1.38243386 

-1.92598991 

-1.97462499 

-2.92349743 

-2.93759426 

-2.15167917 

-2.22825484 

-2.29255366 

-2.36756996 

-2.444138 

-2.5218461 

-2.59716482 

-2.67156189 

-2.74357614 

-2.81235769 

-2.37935869 

-2.94325927 

-3.88499413 

-3.86511716 

-3.12428864 

3.89996478 

3.84855195 

2.97997212 

2.9176184 

2.35391322 

2.78586386 


-99.6433936 

-98.7416369 

-98.3527412 

-99.5384193 

-91.1272886 

-91.2959641 

-91.4899677 

-91.7129879 

-91.9693511 

-92.2648214 

-92.4826918 

-92.9918344 

-93.439864 

-93.9527686 

-94.3427923 

-95.2199546 

-95.9978594 

-96.3831655 

-97.9991336 

-99.077593 

-198.412913 

-191.93694 

-193.66913 

-195.634837 

-197.355598 

-119.351133 

-113.137719 

-116.224333 

-119.613383 

-123.28218 

-127.211277 

-131.353696 

-135.651815 

-149.338842 

-144.445353 

-148.836634 

-153.36923 

-157.19539 

-161.164932 

-164.974791 

-168.636395 

-172.173546 

-175.61834 

-179.988617 

177.614962 

174.210856 

170.739887 

167.167281 

163.465675 

159.613253 
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Figure    9    (Con.t .  ) 


PROBLEM    IDENTIFICATION    -    EX1    W -PLANE   T= . 8 1 

XX    BODE    PLOT    (AMPLITUDE)     XX 


l  I        1      I    I    I  1  I  I 


T 1     I    I    I  I  I  I 


T 1     Mill 


± 


I       ! '    '    tin 


-L 


J I I    1    MM 


J I     I    1    I  1  I 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  = 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  = 
AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +  -60 


.  1  RADIANS/SEC 
108  RADIANS/SEC 
DECIBELS 
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Figure  9  (Cont . ) 


PROBLEM  IDENTIFICATION  -  EX1  W -PLANE  T=  .  8  1 

XX    BODE    PLOT    (PHASE)     XX 


T I 1     I    I    I  I  I  I 


T 1 1     I    I    I  I  I  I 


t 1 — r 


p^rrn 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  98  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  109  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  =  188  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  =  -188  DEGREES 
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FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  EX  1  W' -PLANE  T= .  1 

GAIN—.  8  189414215 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

-2439.71742 
18  1.985871 

1 


NUMERATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 


-121.985871  9 

29  9 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF 

0 

9.24234314 

1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

8  8 

-9.24234314  8 
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Figure    10    (Cont.  )PR0BLEM    IDENTIFICATION    -    EX1    H' -PLANE   T= .  1 
RADIAN  FREQ.       REAL  PART  IMAGINARY  PART     MAGNITUDE  PMSE(RAD)  PWSE(DEG) 


.187226722 

.123284674 

.141747414 

.162975884 

. 187381743 

.215443469 

.247787636 

.284883588 

.327454917 

.376493582 

.43287613 

.497782359 

.572236769 

.657933228 

.756463332 

.369749888 

1.88888881 

1.149757 

1.32194116 

1.51991189 

1.74752841 

2.88923382 

2.31812972 

2.65688781 

3.35385554 

3.51119177 

4.8378173 

4.64158888 

5.33669929 

6.13598734 

7.85488239 

3.1113884 

9.32683358 

18.7226724 

12.3284676 

14.1747418 

16.2975886 

18.7381745 

21.5443472 

24.7787639 

28.4883591 

32.7454921 

37.6493587 

43.2876135 

49.7782365 

57.2236776 

65.7933236 

75.6463342 

86.9749819 

188.888982 


.749981617 
.749869949 


.749772759 

.749699629 

.749682979 

.749475252 

.749386472 

.749883476 

.748788896 

.748399844 

.747886175 

.747288246 

.746313996 

.745135214 

.743582787 

.741548723 

.738858839 

.735343947 

.738749883 

.724766284 

.717889119 

.787812534 

.694228885 

.678841196 

.657796528 

.632868475 

.682754245 

.567282992 

.526356326 

.488862856 

.431916356 

.38117579 

•.338561899 

.281976491 

.237827582 

.196848811 

.161992858 

.  132569482 

•.18827854 

.8885946878 

.8728831383 

.868493815 

.8588148715 

.8433112758 

.9375272875 

.8338884779 

.8296938247 

.8271844597 

.8251332985 


-46.6216872 

-48.5467918 

-35.2628823 

-38.6666494 

-26.6686887 

-23.1989122 

-28.165587 

-17.5334525 

-15.2434888 

-13.2587881 

-11.5164835 

-18.8869112 

-8.6925516 

-7.54777687 

-6.55826855 

-5.63859274 

-4.92181645 

-4.25918654 

-3.6798489 

-3.17268863 

-2.72772751 

-2.33675582 

-1.99239681 

-1.68848289 

-1.41949479 

-1.1813151 

-.978377373 

-.784821224 

-.62833237 

-.478887285 

-.356146355 

-.253986373 

-.178619167 

-.18476951 

-.8546999148 

-.8182671584 

6.98255476E-83 

.8231913619 

.8327938441 

.8375858315 

.8398773822 

.8384219724 

.8364598583 

.3337766111 

.8387785931 

.8276968923 

.8247128768 

.821988753 

.8193382736 

.8169946678 


46.6277179 

48.5537253 

35.2787736 

38.6758137 

26.6792163 

23.2838238 

28.1794298 

17.5494563 

15.2618832 

13.271348 

11.5487753 

18.8348196 

3.72468736 

7.58458359 

6.59251428 

5.72985389 

4.97736475 

4.32279798 

3.75268159 

3.25567824 

2.82237198 

2.44428514 

2.11412199 

1.32555618 

1.57311962 

1.35211887 

1.15851483 

.988939817 

.348554272 

.71181473 

.598388931 

.581859694 

.417619384 

.346767674 

.287233846 

.237738363 

.196968999 

.163644589 

.13656519 

.114616216 

.8968299546 

.8823985323 

.878638582 

.3618164784 

.8531298586 

.8466412816 

.8412985915 

.3369814456 

.8332913884 

.8383397665 


-1.58687981 

-1.58928821 

-1.59285717 

-1.59524863 

-1.59898859 

-1.68318826 

-1.68794548 

-1.61358621 

-1.61989834 

-1.62724569 

-1.63569823 

-1.64539465 

-1.65654519 

-1.66935467 

-1.6848659 

-1.78895514 

-1.72833568 

-1.74256117 

-1.7688284 

-1.79717981 

-1.83849917 

-1.36851612 

-1.91179885 

-1.96889918 

-2.81641584 

-2.97387135 

-2.14878819 

-2.22622829 

-2.31148516 

-2.48429639 

-2.59418884 

-2.61888982 

-2.72872574 

-2.83466462 

-2.94998582 

-3.8646771 

3.18654824 

2.99939685 

2.39989561 

2.38749896 

2.72618143 

2.65644429 

2.59919126 

2.55495334 

2.5238862 

2.58578465 

2.58889615 

2.58593127 

2.52287784 

2.5478319 


-98.9215483 

-91.3595392 

-91.213189 

-91.4395381 

-91.6182882 

-91.3513792 

-92.1285223 

-92.4471289 

-92.8133714 

-93.2343435 

-93.7181381 

-94.274283 

-94.9138817 

-95.6478111 

-96.4899828 

-97.4575356 

-98.5688892 

-99.3414361 

-181.388682 

-182.973889 

-184.879914 

-187.358126 

-189.53754 

-112.351287 

-115.532159 

-119.113597 

-123.111955 

-127.553873 

-132.438391 

-137.756885 

-143.474919 

-149.5426 

-155.386158 

-162.414377 

-169.921752 

-175.593126 

177.991788 

171.352796 

166.186882 

169.85744 

156.198746 

152.2931 

148.922743 

146.388995 

144.688879 

143.578936 

143.245889 

143.579337 

144.584467 

145.93423 
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Figure    10    (Cont . ) 


PROBLEM    IDENTIFICATION    -    EX  1    H'-PLANE   T= . 1 

XX    BODE    PLOT    (AMPLITUDE)     XX 


T 1     MUM 


T 1     11)111 


T 1     I    I    I  II 


I 


J I I     I     It    I  I 


-L 


I       '     '    I    '  '  I  ' 


J_ 


J I      1    t    I  I   I 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 

ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 

MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 

MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  160  RADIANS/SEC 

AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +--49  DECIBELS 
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Figure  10  (Cont . ) 


PROBLEM  IDENTIFICATION  -  EX  1  W'-PLANE  T=.l 

XX  BODE  PLOT  (PHASE)  XX 


T 1 1  I  I  I  I  I  I 


T~T 1  I  I  I  I  II 


T 1 1  I  I  I  I 


J. 


'   I  I 


LLL1 


-L 


'  ■  '  '  ' 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  90  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  109  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  =  188  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  =  -188  DEGREES 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  S-PLANE 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

1 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

0 

19 

1 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 

0  0 

-10  0 

MIN  GAIN  MAX  GAIN 

0  30  0 


Figure  11   s-Plane  Root  Locus  Example 
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Figure  11  (Cont . ) 
1      GAIN  =  9 


GAIN  =  .492857143 


GAIN  =  1.05964234 


GAIN  =  1.71144643 


GAIN  =  2.46102054 


GAIN  =  3.32303076 


GAIN  =  4.31434252 


3      GAIN  =  5.45435104 


GAIN  =  6.76536084 


10     GAIN  =  3.2730221 1 


11     GAIN  =  10  .0063326 


12     GAIN  =  12.0007146 


13     GAIN  =  14.2936739 


ROOTS 

ARE 

REAL  PART 

I  MAG 

-10 

0 

0 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.0495310467 

0 

-9.95046896 

0 

ROOTS 

ARE 

REAL  PART 

I  MA  6 

-.  1071 1 1575 

0 

-9.89288843 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.  174178456 

0 

-9.32582155 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.252476491 

9 

-9.74752351 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.344146776 

0 

-9.65585323 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.451351202 

0 

-9.5481483 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.57395386  1 

0 

-9.42104614 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.729796355 

0 

-9.279  20  365 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-.910  137173 

0 

-9.08936283 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-1.  12739833 

0 

-8.872101 17 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-1.39454782 

0 

-8.60545218 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG 

-1  .72794349 

0 

-8.27205151 

9 
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PART 


PART 


PART 


PART 


PART 


PART 


PART 


PART 


PART 


PART 
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Figure  11  (Cont . ) 

14     GAIN  =  16.938587? 


15     GAIN  =  19.9638332 


16     GAIN  =  23.4583453 


17     GAIN  =  27.4607543 


18     GAIN  =  32.0727246 


19     GAIN  =  37.3764904 


20     GAIN  =  43.4758212 


21     GAIN  =  50 .4900515 


22     GAIN  =  58.5564163 


23     GAIN  =  67.8327359 


24     GAIN  =  78.5005034 


25     GAIN  =  90.7684361 


26     GAIN  =  104.876559 


ROOTS. 

ARE 

REAL 

PART 

I MAG.  PART 

-2.  15932893 

0 

-7.34067107 

0 

ROOTS 

ARE 

REAL 

PART 

I  MAG.  PART 

-2.75568122 

0 

-7.24431878 

0 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-3.75514874 

0 

-6.24435127 

0 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-1  .56867915 

-5 

1  .56367915 
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REAL 
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I MAG.  PART 

-5 

-2.65945945 

-5 

2.65945945 
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ARE 

REAL 

PART 

I MAG.  PART 

-5 

-3.51302365 

-5 

3.51802365 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-4.29835898 

-5 

4.29835898 
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REAL 

PART 

I MAG.  PART 

-5 

-5.84376732 

-5 

5.04876732 
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REAL 

PART 

I MAG.  PART 

-5 

-5.79279803 

-5 

5.79279883 

ROOTS 

ARE 
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I MAG.  PART 

-5 

-6.54467233 

-5 

6.54467233 
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I MAG.  PART 

-5 

-7.31448383 

-5 
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I MAG.  PART 

-5 
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REAL 

PART 

I MAG.  PART 

-5 

-3.93736367 

-5 
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Figure  11  (Cont .  ) 

27     GAIN  =  121.  1889 


28     GAIN  =  139.758392 


29     GAIN  =  161.215582 


38     GAIN  =  185.898777 


31     GAIN  =  214.267251 


32     GAIN  =  246.988 195 


33     GAIN  =  284.423882 


ROOTS  ARE 


REAL 

-5 

-5 

REAL 
-5 

-5 

REAL 

-5 

-5 

REAL 
-5 


REAL 

-5 

-5 

REAL 

-5 

-5 

REAL 

-5 

-5 


PART 


ROOTS  ARE 


PART 


ROOTS  ARE 


PART 


ROOTS  ARE 


PART 


ROOTS  ARE 


PART 


ROOTS  ARE 


PART 


ROOTS  ARE 


PART 


I MAG.  PART 

-9.88318663 

9.38318663 

I MAG.  PART 
-18.7125577 
18  .7125577 

I MAG.  PART 
-11.6711431 
1  1  .671 1431 

I MAG.  PART 

-12.6342723 

12.6342728 

I MAG.  PART 

-13.7574435 

13.7574435 

I MAG.  PART 

-14.3963148 

14.8963148 

I MAG.  PART 
-16.  18  67713 
16.  1367713 


Figure  11 
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Figure  11  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EX  1  S-PLANE 
XX  ROOT  LOCUS  PLOT  XX 


I  1  I  I  11  I  I  I  I  I  I  I  I  i  I  I  I  I  1  I 


+ 
+ 
* 
♦ 

+ 


I  I  I  I  I  I  I  I  II 


i  i  i  t  i  i  i  i  i  i  i  <»  i  i  i  I  i  I  i  » 


'«'''■'''' 


ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHON  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -22  TO  1  1 

ORDINATE,  -17  TO  16 


Figure  11 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  W -PLANE  T=.88  1 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

-2.485638E+8? 

128831? 

1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 


-128281?  8 

2888  8 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 


9 

?, 
1 

,???? 

Eh 

* 

J-LOOP  POLES 
REAL  PART 

8 
-?.???? 

MIN  GAIN 
8 

IMAGINARY  PART 


9 

8 

MAX 

GAIN 

-1  . 

2546E- 

-07 

0*SXX***S******X*******XX*****XS**********X****** 


Figure  12   w1 -Plane  Root  Locus  Example,  T=.001 
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Figure  12  (Cont . ) 
1      GAIN  =  8 


GAIN  =  -2.961  12857E-10 


GAIN  =  -4.43142643E-10 


GAIN  =  -7. 15726896E-10 


GAIN  =  -1.02919879E-99 


GAIN  =  -1  . 33969  146E-99 


GAIN  =  -1.30425804E-09 


8 


10 


11 


12 


13 


GAIN  =  -2.2310096E-09 


GAIN  =  -2.3292739E-99 


GAIN  =  -3.45977784E-09 


GAIN  =  -4.  13435737E-09 


GAIN  =  -5.0  1369884E-09 


GAIN  =  -5.97761652E-09 


ROOTS  ARE 
REAL  PART  I MAG.  PART 

-9.9999  0 

0  0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.0498333602       9 
-9.94981915        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.107772379        0 
-9.89159549        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.175266637        9 
-9.82377392        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.254077448        0 
-9.74458669        8 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.346367783        0 
-9.65186347        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.454329223        9 
-9.54299422        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.53236724         9 
-9.4  1429371        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

- . 734883239        9 
-9.26161935        9 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-.916725984        9 
-9.97901949        0 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-1.13646232        0 
-8.35841248        0 

ROOTS  ARE 
REAL  PART  IMAG.  PART 

-1.40580821        0 
-3.53806531        0 

ROOTS  ARE 

REAL    PART  IMAG.     PART 

-1.74311005  0 

-8.24961193  0 

Figure    12 
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Figure  12  (Cont 
14     GAIN  =  -7 


.) 
.03037185E-09 


15     GAIN  =  -S.34854049E-09 


16     GAIN  =  -9.30693442E-09 


17     GAIN  =  -1.  14840874E-03 


13     GAIN  =  -1 .34123134E-98 


19     GAIN  =  -1 .56308483E-08 


20     GAIN  =  -1.31815884E-03 


21     GAIN  =  -2.  1  1  149395E-03 


22     GAIN  =  -2.44882933E-08 


23     GAIN  =  -2.3367650  1E-08 


24     GAIN  =  -3.23289105E-08 


25     GAIN  =  -3.79593599E-08 


26     GAIN  =  -4.33593767E-03 


ROOTS  ARE 
REAL  PART 
-2.  18070408 
-7.81069376 

ROOTS  ARE 
REAL  PART 
-2.7390739 
-7.2008011 

ROOTS  ARE 
REAL  PART 
-3.83271797 
-6. 15540578 

ROOTS  ARE 
REAL  PART 
-4.9930549 
-4.9930549 

ROOTS  ARE 
REAL  PART 
-4.99139689 
-4.99139639 

ROOTS  ARE 
REAL  PART 
-4.99056517 
-4.99056517 

ROOTS  ARE 
REAL  PART 
-4.98903369 
-4.93903369 

ROOTS  ARE 
REAL  PART 
-4.9872725 
-4.9872725 

ROOTS  ARE 
REAL  PART 
-4.93524712 
-4.93524712 

ROOTS  ARE 
REAL  PART 
-4.9829  1794 
-4.98291794 

ROOTS  ARE 
REAL  PART 
-4.98023937 
-4.98023937 

ROOTS  AF 
REAL  PART 
-4.97715903 
-4.97715903 

ROOTS  AF 
REAL  PART 
-4.97361663 
-4.97361663 
Figure  12 
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I MAG.  PART 

0 

0 

I MAG.  PART 

0 

0 

I MAG.  PART 
0 

0 

I MAG.  PART 
-1  .64  193793 
1  .64  193793 

I MAG.  PART 
-2.71059762 
2.7  1059762 

I MAG.  PART 
-3.56320  396 
3.56320396 

I MAG.  PART 
-4.34  14  1261 
4.34  14  1261 

I MAG.  PART 
-5.09  13667 
5.0913667 

I MAG.  PART 

-5.33586167 

5.83586167 

I MAG.  PART 
-6.583341 16 
6.583841 16 

I MAG.  PART 
-7.360 14213 
7.360 14218 

I MAG.  PART 
-8.  15747327 
8.  15747327 

I MAG.  PART 
-8.9873759 
3.9873759 


Figure  12  (Cont . ) 

27     GAIN  =  -5.06443961E-08 


28     GAIN  =  -5.34471684E-08 


2?     GAIN  =  -6.74203565E-08 


30     GAIN  =  -7.77395223E-08 


31     GAIN  =  -8. 9406564  1E-08 


32     GAIN  =  -1  .03253662E-0 


33     GAIN  =  -1. 13947324E-07 


ROOTS  ARE 
REAL  PART 
-4.96954288 
-4.96954288 

ROOTS  ARE 
REAL  PART 
-4.96485806 
-4.96485806 

ROOTS  ARE 
REAL  PART 
-4.95947052 
-4.95947052 

ROOTS  ARE 
REAL  PART 
-4.95327434 
-4.95327484 

ROOTS  ARE 
REAL  PART 
-4.94614981 
-4.94614981 

ROOTS  ARE 
REAL  PART 
-4.93795603 
-4.9379560  3 

ROOTS  ARE 
REAL  PART 
-4.92353319 
-4.92853319 


I MAG.  PART 

-9.35574623 

9.85574623 

I MAG.  PART 

-10.7631438 

10.7681433 

I MAG.  PART 
-1  1.7299837 
1 1 .7299387 

I MAG.  PART 
-12.746695 
12.746695 

I MAG.  PART 

-13.8237636 

13.3237636 

I MAG.  PART 

-14.9663317 

14.^668817 

I MAG.  PART 
-16. 1319344 
16.  1819344 


Figure  12 
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Figure  12  (Cont . ) 


PROBLEM  IDENTIFICATION  -  EX  1  H' -PLANE  T= 
XX    ROOT  LOCUS  PLOT  XX 


801 


I  I  I  I  I  I  I  I  I  I  I  !  I  I  I  I  I  I  I  I  I 


+ 

+ 
♦ 

4 

■f 
+ 


I  I  I  I  I  I  I  II  I 


II  I  I 


♦  I  IIM 


I  I  1  I  I  I  I  I  I  I  I  I  I  I  I  1  I  I  I  I  1 


I  I  I  I  I  I  1  I  1  1 


ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -22  TO  11 

ORDINATE,  -17  TO  16 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  W -PLANE  T= . 8 1 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W' 

-2480409.84 
1 1802.8042 
1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

200  8 

-12802.0042  8 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W' 


0 

9.99167498 
1 

OPEN-LOOP  POLES 
REAL  PART 

IMAGINARY 

0 

-9.99167498 

0 

8 

MIN  GAIN 

0 

MAX  GAIN 
-1  .25E-84 

************************************************** 


Figure  13   w1 -Plane  Root  Locus  Example,  T=.01 
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Figure  13  (Cont . ) 
1      GAIN  =  9 


GAIN  =  -2.95357143E-97 


GAIN  =  -4.41517357E-97 


GAIN   =    -7.  13182678E-07 


GAIN  =  -1  .92542522E-06 


GAIN  =  -1. 384596 15E-06 


GAIN  =  -1.79764271E-06 


3      GAIN  =  -2.27264626E-06 


GAIN  =  -2.81390035E-06 


10     GAIN  =  -3.44709254E-06 


11     GAIN  =  -4.  16951357E-06 


12     GAIN  =  -5.90029774E-06 


13     GAIN  =  -5.95569955E-06 


ROOTS  ARE 
REAL  PART  I MAG.  PART 

-9.99167498        9 

0  0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.0495932003       0 
-9.93966022        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.107273051        0 
-9.87919055        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.174510  775        0 
-9.30375529        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.253064237        0 
-9.72651892        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.345120662        0 
-9.63022713        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.453403097        0 
-9.51707403        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.581368364        0 
-9.38359749        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.73350241         9 
-9.22493198        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.915796521        0 
-9.03523017        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1.1365675         0 
-8.80594033        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1.40799799        0 
-8.52471313        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1.74939626        9 

-8.17294862        9 

Figure  13 
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Figure  13  (Cont 
14     GAIN  =  -7 


15     GAIN  =  -8 


16     GAIN  =  -9 


17     GAIN  =  -1 


18     GAIN  =  -1 


1?     GAIN  =  -1. 


20     GAIN  =  -1 . 


21     GAIN  =  -2. 


22     GAIN  =  -2. 


23     GAIN  =  -2 


24     GAIN  =  -3 


25     GAIN  =  -3 


26     GAIN  =  -4 


.) 

.85441 162E-06 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-2.  1954399 

9 

-7.71304879 

0 

.31793051E-06 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-2.82447443 

0 

-7.0691 1459 

0 

.77097722E-06 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-3.97303464 

0 

-5.90341974 

0 

.  1441981E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.92837473 

-1  .73235163 

-4.92837473 

1  .73235163 

.33636352E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.91704436 

-2.8109046 

-4.91704436 

2.3109046 

■55735377E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.90401439 

-3.65156925 

-4.9040 1439 

3.65156925 

.81149255E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.33902985 

-4.42507143 

-4.88902985 

4.42507143 

.  10375214E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.37179752 

-5.  17350  32 

-4.37179752 

5.  17350  32 

.43985068E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.35198023 

-5.91326942 

-4.85198023 

5.91826942 

.32636399E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.32919018 

-6.67269846 

-4.829190  18 

6.67269846 

.  2708543E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.30293141 

-7.4463113 

-4.80298141 

7.44631 13 

.7820  1816E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.77234103 

-8.24663677 

-4.77284103 

3.24663677 

.  3698566E-05 

ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-4.73817921 

-9.030 10553 

-4.73817921 
Figure  13 
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Figure  13  (Cont . ) 

27     GAIN  =  -5.94587931E-95 


28     GAIN  =  -5.32323715E-95 


29     GAIN  =  -6.71731593E-95 


39     GAIN  =  -7.74544993E-95 


31     GAIN  =  -8.9273021E-05 


32     GAIN  =  -1  .02875081E-04 


33     GAIN  =  -1  .  1851  179  1E-94 


ROOTS  ARE 
REAL  PART 
-4.69331762 
-4.69331762 

ROOTS  ARE 
REAL  PART 
-4.65247613 
-4.65247613 

ROOTS  ARE 
REAL  PART 
-4.59975752 
-4.59975752 

ROOTS  ARE 
REAL  PART 
-4.53912996 
-4.53912996 

ROOTS  ARE 
REAL  PART 
-4.46940673 
-4.46940673 

ROOTS  ARE 
REAL  PART 
-4.38922296 
-4.38922296 

ROOTS  ARE 
REAL  PART 
-4.29700894 
-4.29790894 


I MAG.  PART 
-9.95253982 
9 . 95253932 

I MAG.  PART 
-19 .3*94452 
19 .3694452 

I MAG.  PART 
-11.3361975 

1  1  .336  1975 

I  MAG.  PrtRT 

-12.3581696 

12.3581696 

I MAG.  PART 
-13.940324 
13.940824 

I MAG.  PART 
-15.9  397333 
15.9S97333 

I MAG.  PART 

-16.3198365 

16.3198865 


Figure  13 
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Figure  13  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  EX1  W -PLANE  T= .  8  1 
XX  ROOT  LOCUS  PLOT  XX 


I  I  I  I  I  I  I  1  I  I  M  I  I  I  1  I  I  M  I  I  I  I  I  I  I!  I  I  I  I 


i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  r:  i  i  i 


I  I  I  I  1  I  I  i  1  I 


ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -22  TO  1 1 

ORDINATE,  -17  TO  16 


Figure  13 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  W'-PLANE  T= .  1 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

-2439.71742 
18  1  .985871 

1 


OPEN-LOOP  2EROS 

REAL  PART  IMAGINARY  PART 


-121.985871  8 

28  8 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W' 

8 

9.24234314 

1 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 

8  8 

-9.24234314  8 

MIN  GAIN  MAX  GAIN 

8  -.1  13643529 


Figure  14   w' -Plane  Root  Locus  Example,  T=.l 
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Figure  14  (Cont .  ) 


1 


5 


3 


10 


11 


12 


13 


GAIN  =  9 


GAIN  =  -1 .8670329SE-04 


GAIN  -  -4.0 142284E-04 


GAIN  =  -6.48344563E-04 


GAIN  =  -9.32304545E-04 


GAIN  =  -1.25835852E-03 


GAIN  =  -1  .6343956E-03 


GAIN  =  -2.06626324E-93 


GAIN  =  -2.56291 102E-03 


GAIN  =  -3.  13405597E-03 


GAIN  =  -3.79937267E-9  3 


GAIN  =  -4.54621 136E-9 3 


GAIN  =  -5.4  1435194E-03 


ROOTS  ARE 
REAL  PART  I MAG.  PART 

-9.24234314        9 
0  9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.0496547373       9 
-9.17536918        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.197695744        9 
-9.99749393        0 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.175741467        0 
-9.9964329  1        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.255897277        9 
-8.89993997        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.359442074        9 
-3.77590293        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.462932327        9 
-8.62753256        9 

ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-.597624754        9 
-3.45268997        9 

ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-.76945195?        9 
-3.24358792        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-.959374979        9 
-7.99039184        9 

ROOTS  ARE 
REAL  PART  I  MAG.  PART 

-1.29372442        9 
-7.63979  196        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1.52831341        9 
-7.29947243        9 

ROOTS  ARE 
REAL  PART  I MAG.  PART 

-1.95977976        9 
-6.77764531        9 
Figure  14 
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Figure  14 
14     GAIN 


(Cont .  ) 
=  -6.4137S803E-03 


15     GAIN  =  -7.56256453E-03 


16     GAIN  =  -8.3836575E-03 


17     GAIN  =  -.9  194929144 


13     GAIN  =  -.0  12150059? 


1?     GAIN  =  -.0  141592772 


20     GAIN  =  -.0  16469377 


21     GAIN  =  -.0 191270669 


22     GAIN  =  -.0221828352 


23     GAIN  =  -.0256969638 


24     GAIN  =  -.9297332224 


25     GAIN  =  -.034385664 


26     GAIN  =  -.039730222 


ROOTS  ARE 
REAL  PART 
-2.61028387 
-6.03338216 

ROOTS  ARE 
REAL  PART 
-4.26780979 
-4.26789979 

ROOTS  ARE 
REAL  PART 
-4.20552327 
-4.29552327 

ROOTS  ARE 
REAL  PART 
-4.  13369394 
-4.  13369394 

ROOTS  ARE 
REAL  PART 
-4.05082203 
-4.05882203 

ROOTS  ARE 
REAL  PART 
-3.95515054 
-3.95515054 

ROOTS  ARE 
REAL  PART 
-3.34464502 
-3.84464502 

ROOTS  ARE 
REAL  PART 
-3.71692006 
-3.71692006 

ROOTS  ARE 
REAL  PART 
-3.56917313 
-3.569  17313 

ROOTS  ARE 
REAL  PART 
-3.39812932 
-3.39812932 

ROOTS  ARE 
REAL  PART 
-3.  1998914 
-3.  1993914 

ROOTS  ARE 
REAL  PART 
-2.96936645 
-2.96986645 

ROOTS  ARE 
REAL  PART 
-2.70258524 
-2.70258524 
Figure  14 

127 


I MAG.  PART 

0 

9 

I MAG.  PART 
-.6  13935262 

I MAG.  PART 

-2.944S5929 

2.94435929 

I MAG.  PART 

-2.92566367 

2.92566367 
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C.   ROOT  LOCUS  TEMPLATES  IN  w » -PLANE 

The  s-plane  consists  of  an  abscissa  representing  the 
real  portion  of  the  complex  variable  s  and  an  ordinate 
representing  the  imaginary  part  where 

s=  0)    +  JU) 

The  w' -plane  approaches  the  s-plane  as  the  sampling 
period  T  approaches  zero.  (See  [Ref.  11]  for  a  proof.) 
The  w1 -plane,  like  the  s-plane,  consists  of  an  abscissa 
representing  the  real  part  and  an  ordinate  representing 
the  imaginary  part  where 

The  s  and  w*  planes  are  related  as  shown  below. 

u1  =  2/T  tanh  (  (J  T/2) 

P%    =  2/T  tan(  (jj '  T/2) 

(J    =   2/T  tanh-1  (u'T/2)     where  -1  <  u'T/2  <  1 

6Jj  =  2/5  tan"1  (  y'T/2)    where  -1  <  J/*T/2  <  1 

In  the  s-plane,  (J   represents  damping.   Therefore  a  line 
of  constant  damping  in  the  s-plane  is  a  vertical  line 
parallel  to  the  imaginary  axis  and  passing  through  the  proper 
value  of  (J    on  the  real  axis.   Since  the  real  axis  in  the 
w' -plane  represents  u'  which  is  not  the  damping  but  only 
related  to  C    as  shown  above,  it  is  not  possible  to  easily 
determine  the  damping  of  a  particular  root  by  simply  ob- 
serving the  position  of  the  root  in  the  plane.   This  is  also 
true  of  the  damped  natural  frequency,  the  natural  frequency, 
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and  the  damping  ratio.   Because  of  this  it  is  helpful  to 
have  templates  of  constant  parameters  in  the  w' -plane  to 
make  interpretation  of  the  characteristics  of  a  root 
by  its  location  easier.   The  following  sections  give  some 
insight  into  the  nature  of  these  templates  in  the  w1 -plane. 
1 .   Constant  Damping 

The  constant  damping  templates  are  shown  in  Figure  15 
There  are  three  templates,  one  of  each  of  three  values  of 
sampling  period,  .001  seconds,  .01  seconds,  and  .1  seconds. 
Each  template  shows  lines  of  constant  damping  for  the  values 
25,  50,  75,  and  100.   These  templates  are  created  by  using 
the  above  relationship  between  (J  and  u'  to  find  the  constant 
value  of  u'  that  corresponds  to  the  constant  value  of  (J    and 
then  plotting  the  constant  u1  value.   From  the  templates  it 
is  seen  that  for  a  period  of  .001  there  is  practically  no 
difference  from  the  s-plane.   For  a  period  of  .01  seconds, 
some  distortion  is  noticed  as  the  values  of  damping  approach 
the  reciprocal  of  the  value  of  the  period.   This  distortion 
is  in  the  form  of  u1  becoming  smaller  than  the  value  of 
damping.   And  for  a  period  of  .1  seconds  gross  distortion 
is  seen.   Figure  15  shows  an  extra  template  for  a  period 
of  .1  seconds  but  for  values  of  damping  of  1,  2,  3,  and  4. 
In  this  case  only  little  distortion  is  present  because  the 
values  of  damping  are  small  compared  to  the  reciprocal  of 
the  period. 
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2 .  Constant  Damped  Natural  Frequency 

Figure  16  shows  the  templates  for  constant  damped 
natural  frequency  for  values  of  damped  natural  frequency  of 
25,  50,  75,  and  100.   The  results  here  are  similar  to  the 
constant  damping  case.   As  the  period  increases,  the 
distortion  between  the  value  of  OJj  and  "the  value  of  V    in- 
creases.  In  this  case  the  value  of  V1   becomes  greater  than 
the  value  of  damped  natural  frequency  as  it  distorts. 
In  the  case  where  the  period  is  .1  seconds  the  distortion 
is  so  great  that  the  template  is  useless  and  negative 
values  of  Pl   are  produced  by  the  tangent  function  in  the 
relationship. 

3 .  Constant  Natural  Frequency 

Templates  for  constant  natural  frequency  in  the 
w1 -plane  are  shown  in  Figure  17.   In  the  s-plane  constant 
natural  frequency  QJn    plots  as  a  circle  with  its  center  at 
the  origin  where 

6Jh  =  sqrt(  C2  +  CJoJ2  ) 
To  plot  constant  natural  frequency  in  the  w1 -plane,  substitute 

u'  =  2/T  tanh(  (J   T/2) 

V  =  2/T  tan(  CO4   T/2) 

into  the  above  equation  and  solve  for  V1   to  get 

V  =   2/T  tan 


T/2  {(J„  -  [2/T  tanh  -1  (u'T/2)]2}* 
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For  a  constant  value  of  natural  frequency  and  period  this 
equation  is  used  to  plot  u1  vs  pl   .   In  the  w'-plane  for  a 
period  of  .001  seconds  and  a  value  of  natural  frequency  of  25 
it  plots  very  close  to  a  perfect  circle.   As  the  values  of 
natural  frequency  increase  they  form  concentric  ellipses 
which  are  elongated  along  the  imaginary  axis  and  contracted 
along  the  real  axis.   This  effect  becomes  more  dramatic 
as  the  period  increases.   An  extra  template  for  a  period  of 
.05  seconds  is  added  to  better  show  the  trend  since  the 
change  from  .01  seconds  to  .1  seconds  is  so  great. 
4.   Constant  Damping  Ratio 

The  damping  ratio,  £     ,  can  be  defined  as 


£=  \CT  I  Un   I  =  cos# 


where     represents  the  absolute  value  and  Q    is  the  angle 
formed  by  the  negative  real  axis  and  a  line  joining  the 
origin  and  the  root  in  a  root  locus  plot.   In  the  s-plane  a 
plot  of  constant  damping  ratio  is  a  line  radially  out  from 
the  origin  forming  the  correct  angle  $  with  the  negative 
real  axis.   To  plot  constant  damping  ratios  in  the  w'-plane, 
substitute  the  proper  relations  into  the  above  equation  to 
get 

tanh"1  (u'T/2) 
^    {[tanh"1  (u'T/2)]2  +  [tan-1  (^'T/2)]2}* 
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Solving  for  u'  gives 


u1  =  2/T  tanh 


tan-1  (  VK    T/2) £ 
sqrtd-  £2  ) 


This  equation  is  used  to  plot  constant  damping  ratio  for 
a  given  sampling  period.   Figure  18  shows  templates  for 
constant  damping  ratio  for  values  of  damping  ratio  of  .1, 
.5,  .707,  and  .9.   From  these  templates  it  can  be  seen  that 
there  is  a  region  near  the  origin  that  can  be  interpreted 
the  same  as  the  s-plane  for  a  close  approximation  and  that 
the  size  of  this  region  depends  on  the  value  of  the  sampling 
period.   Beyond  this  region  a  template  makes  the  interpreta- 
tion easier. 
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Figure  15   Constant  Damping  Templates 
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Figure  15  (Cont . ) 
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Figure  15  (Cont . ) 
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Figure  15 
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Figure  16  (Cont . ) 


CONSTANT  DAMPED  NATURAL  FREQ  N' -PLANE  T= . 0  1 
XX  ROOT  LOCUS  PLOT  XX 


5  ii  i  1 1  1 1  ii  i  rrn  n  1 1  r 

j  i  i  r  1 1  ii  ii  i 

— 

- 

■ 

- 

- 

- 

- 

- 

Xi.1.1  I1IILIIII.1MIIII 

«i i  i  (  i  i  i  i  i  i  i  »  i  r 

ABSCISSA  ->  REAL  < U)  AXIS 
ORDINATE  ->  IMAGINARY  (NU)  AXIS 
TIC  MARKS  SHON  INTERNALS  OF  5 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -1O0  TO  100 
ORDINATE,  -32  TO  163 


Figure  16 
140 


Figure  16  (Cont . ) 
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Figure  17   Constant  Natural  Frequency  Templates 
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Figure  17  (Cont . ) 
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Figure  17  (Cont . ) 
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Figure  17  (Cont . ) 
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Figure  18   Constant  Damping  Ratio  Templates 
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Figure  18  (Cont.  ) 
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J  n  I  I  I  I  II  l:  I  I  I  I  I  II  M  I  I  I  I  I  I  II  .4  I 


i  i  i  i  i  i  i  i  i 


\     % 


t  1  I  1  I  1  I  1  I  1  1  1  I  I  I  t  1  i    I  I  I  I  I  I  t  I  I  f  i- 


I  I  I  I  I  I  I  I  I 


ABSCISSA  ->  REAL  < U)  AXIS 
ORDINATE  ->  IMAGINARY  (NU)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  10 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -369  TO  100 

ORDINATE,  0  TO  400 


Figure  18 
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Figure  18  (Cont .  ) 


CONSTANT  DAMPING  RATIO  W -PLANE  T= . 1 
XX  ROOT  LOCUS  PLOT  XX 


J  11  I  I  I  I  I  I  l:ii  I  I  i  I  I  !  }    !  I  1  I  «  I  I  1  )   I 


I  I  I  I  I  I  I  I  I 


»  I  I  »»  I  I  I  I  I  t  I  I  I  I  1  1  I  |  M  I  (  I  I  I  I  <  i- 


I  t  I  I  1  I  I  I  I 


ABSCISSA  ->  REAL  (U)  AXIS 
ORDINATE  ->  IMAGINARY  <NU)  AXIS 
TIC  MARKS  SHON  INTERNALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -36  TO  18 

ORDINATE,  0  TO  40 


Figure  18 
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D.   COMPENSATION  EXAMPLE 

To  further  demonstrate  the  use  of  the  RTLOC  and  FRESP 
programs  for  analysis  and  design  in  the  w' -plane,  an  example 
from  [Ref.  6]  is  duplicated  here. 

An  angular  position  servo  is  described  below. 


Md  (friction) 


Mf3%->' 


K  D 


t>ZOH 


P 


>K. 


-^Kh?&— > 


Js' 


> 


where 


K   =  1.9  (V/rad) 
P 

KT  =  0.385  (Nm/A) 


K,  =  KA/V) 

J  =  3.22  x  103   (Nm.  Sec2) 


Neglecting  friction  (Md)  gives 


ft— ^5^0-^-^V 


KIKT    l-e"sT 


■» 


P  ^ 


where 


1-e 


-sT 


ZOH  = 


149 


-sT 
KTKT   1-e 

G1(S)  =  ~-       —^ 


Converting  to  the  z  domain  gives 

(z+1) 


G,  (z)  =  K0T2 


2    r       in2 
(z-1) 


where 


K„  =  — — -  =59.8 
Z  2J 


Converting  this  to  the  w  domain  gives 


K9T2    1-w 
2      w 


where 


w  =  w'T/2 


The  open  loop  transfer  function  is 


G(w)  =  K  K.G.  (w) 
p  1  1 


where 


Kx  =  3.83 
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from  the  calculation  of  D  C  gain.   This  gives 


9    1-w 
G(w)    =   217.58    T      — j- 

w 


For 


T   =   0.01 

1-w 
G(w)    =   0.0218   ~y  (5) 

w 


Therefore  for  the  uncompensated  system  equation  5  above  is 
entered  into  the  FRESP  program  producing  the  output  shown 
in  Figure  19.   From  these  Bode  plots  it  is  seen  that  the 
uncompensated  system  is  unstable. 

From  analog  design  techniques  the  final  compensation 
network  [Ref.  7]  is 

(w/0.1  +  1) 

K-D(w)  =  3.83  P 

(w/2.15  +  1) 

and  the  open  loop  transfer  function  becomes 

w+0.1     w-1 
G(w)  =  K  


(w+2.15)2    w2 


For 


K  =  -1.017 
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this  transfer  function  was  entered  into  the  FRESP  program. 
The  resulting  Bode  plots  in  Figure  20  show  that  the  compen- 
sated system  is  stable  with  a  gain  margin  of  10.5  dB  and 
phase  margin  of  41  degrees. 

Figure  21  shows  the  output  from  the  RTLOC  program  for  the 
compensated  system.   This  figure  includes  the  normal  root 
locus  in  addition  to  an  expanded  portion  of  the  root  locus 
to  obtain  more  detail.   These  plots  are  of  w  and  not  w'  in 
this  case.   Values  of  w1  can  be  found  by  the  simple  relation 

w1  =  2w/T 
where 

T  =  .01  seconds 
for  this  problem. 

On  the  expanded  root  locus  plot  in  Figure  21,  a  root 
associated  with  a  certain  gain  may  be  identified  and 
templates  similar  to  the  ones  discussed  in  section  V.C.  may 
be  used  to  easily  obtain  the  characteristics  associated  with 
the  root. 
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FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

X**XXX**** ******* *£***S*SXXX*******£*£S*XXX**SX£X: 

GAIN=-.822 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

-1 
1 


NUMERATOR  ROOTS  ARE 

'  REAL  PART  IMAGINARY  PART 

1  9 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

9 
@ 

1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PhRT 

9  9 

9  9 


Figure  19   Uncompensated  System  Frequency  Response 
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PROBL EM  I D  EMT I F I  CAT  I  ON 

-  KATZ 

EXAMPLE 

RADIAN  FREQ. 

REAL  PART 

IMAGINARY  PART 

MAGNITUDE 

PMSE(RAD) 

PHASE(DEG) 

.194761575 

-2.99455394 

.21889944 

2.31552385 

3.33721183 

174.319434 

.1149757 

-1.44421932 

.19134478 

1.4751832 

3.92711966 

173.441242 

.126135688 

-1.33144411 

.  174344238 

1.39242847 

3.91687943 

172.333163 

.133488637 

-1.14783132 

.158857798 

1.15382958 

3.99397933 

172.115482 

.151991199 

-.952327477 

.  144745389 

.943244472 

2.99875693 

.164819954 

-.798439891 

.131334535 

.881543527 

2.97639452 

179.529743 

. 183973329 

-.656482338 

.129179998 

.447311733 

2.968524 

169.425591 

.299923391 

-.544956794 

.199494513 

.555847942 

2.94339951 

163.o39273 

.229513975 

-.452432784 

.9997473267 

.463382123 

2.92455396 

167.544687 

.242912328 

-.375417579 

.8989842724 

.386441842 

2.99414539 

.  -J  .  -  7~.  -  -  J 

.24560878 

-.311344312 

.3323235873 

.322454841 

2.33197322 

165. 125243 

.291595393 

-.253393427 

.8754793239 

.249474375 

2.35794731 

163.743377 

.319926716 

-.214942188 

.3437457432 

.225474241 

2. S3 195624 

162.259393 

.351119176 

-.17344373 

.8424547834 

.139129143 

2.39392119 

169.652903 

.385352342 

-.148151312 

.3579985323 

.153778714 

2.77377459 

153.925743 

.422924291 

-.122997322 

.3523137447 

.133545559 

2.74143151 

157.975374 

.444158337 

-.182114953 

.3473975428 

.112573323 

2.79732683 

155.191263 

.539413394 

-.3847774274 

.3431348939 

.8951433593 

2.67944244 

153.335126 

.559981923 

-.8783838757 

.8393582893 

.3884378584 

2.63138424 

159.791329 

.413599733 

-.8584339391 

.835354513 

.8435578413 

2.5912399 

143.467163 

.473415972 

-.8485123753 

.8324493315 

.3584874545 

2.54393264 

146.843126 

.739972211 

-.3492742413 

.3297473587 

.3599824412 

2.59512214 

143.532977 

.31113984 

-.8334389431 

.827122428 

.9439559773 

2.46913142 

149.952479 

.893215995 

-.9277688598 

.3247131234 

.9371472377 

2.41421 

123.324993 

.977999949 

-.323947544 

.3225174322 

.8322214493 

0  Oi'^OnO"7 
CtCOtOiX/ 

135.444294 

1.97224724 

-.9191344774 

.9285172734 

.9238559462 

2.32133514 

133.382755 

1.17481197 

-.8158357794 

.3134945753 

.3245325322 

2.27514727 

lco .--.--. 

1.29154948 

-.3131334531 

.3179333883 

.32154277 

2.22944425 

127.749364 

1.41747418 

-.9199494515 

.3155285449 

.9189941431 

2.18519983 

125.292259 

1.55547417 

-9.39942792E-33 

.3141417497 

.3148114621 

2.14299359 

122.733252 

1.79735247 

-7. 5479322 1E-93 

.3123854454 

.3149329282 

2.18963635 

129.35764 

1.37381745 

-6.24547372E-93 

.3117497281 

.31338383 

2.36131571 

113.337544 

2.85451234 

-5.29137454E-93 

.9196977233 

.311395411 

2.32339131 

115.931352 

2.25781976 

-4.31343929E-93 

9.74734794E-93 

.8134412495 

1.93736531 

113.394342 

2.4779764 

-3.53545172E-93 

3.83143732E-93 

9.57735993E-93 

1.95448931 

111.934957 

2.71853329 

-2. 9747959 1E-93 

3.39243536E-93 

3.42254583E-93 

1.92327327 

113.195431 

2.93344729 

-2.47131279E-83 

7.3735257E-93 

7.77444891E-93 

1.39419326 

183.529144 

3.27454922 

-2.95172729E-93 

4.71848199E-93 

7.92473344E-83 

1.36713643 

194.931955 

3.59281373 

-1.78333984E-83 

4.12143957E-83 

4.35429953E-93 

1.34213798 

135.549532 

3.94428413 

-1.4H17483E-93 

5.57789 173E-93 

5.75428244E-93 

1.31918968 

194.224329 

4.32374137 

-1.17497378E-83 

5.38223524E-93 

5.21413573E-93 

1.79732442 

193.387915 

4.75981924 

-9.74734482E-84 

4.43873992E-83 

4.73224353E-93 

1.77325322 

191.334727 

5.2148884 

-8. 992435 14E-84 

4.21948249E-33 

4. 2943845 1E-83 

1.76928437 

189.357914 

5.72234778 

-4.71843134E-94 

3.3445424E-93 

3.99282442E-33 

1.74389221 

99.9125428 

4.23929158 

-5.57739 164E-94 

3.58332287E-83 

3.5471512E-93 

1.72349999 

99.347197 

4.39261226 

-4.6387339 1E-84 

3.19132324E-93 

3.22524372E-83 

1.71437394 

93.2559754 

7.56443344 

-3. 3445623 1E-94 

2.99827844E-83 

2.93357 183E-83 

1.79222342 

97.5395339 

3.39217537 

-3.19132317E-94 

2.44998773E-93 

2.46996132E-83 

1.49944924 

94.3432463 

9.11142777 

-2.64998747E-84 

2.41449723E-83 

2.42899585E-93 

1.63919374 

96.2431754 

19.9899992 

-2.19999992E-34 

2.19999995E-93 
Fig 
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Figure  19  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

XX  BODE  PLOT  (AMPLITUDE)  XX 


T 1 1     i  I  I  1 


i  i  i  i  r 


! 


I 


_L 


_L 


J 1111 


J I   I   1 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 
MINIMUM  FREQUENCY  SHOWN  ON  hBSCISSA  = 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  = 


.  1  RADIANS/SEC 
10  RADIANS/SEC 


AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-69  DECIBELS 


Figure  19 
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Figure  19  (Cont. ) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

XX  BODE  PLOT  (PHASE)  XX 


1 1 1 — I  I  I  I 


T 1 1 1 — I  I  I 


_L 


J I I   I  I  1 


_L 


-L 


I  \     I  I 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHON  MULTIPLES  OF  99  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  18  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  =  180  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  =  90  DEGREES 


Figure  19 
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F REG U EIM C Y  R  E  S  P  ON  S E 

PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

GAIN=-1 .9 17 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  : 

-.  1 
-.? 

1 


NUMERATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

-.1  9 

1  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

0 

0 

4  .  6225 

4.3 

1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

-2. 15  0 

-2.15  0 

0  0 

0  0 


Figure  20   Compensated  System  Frequency  Response 
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Figure  20  (Cont .  ) 


P  R  0  B  L  EM    I  D  ENT I  F I  CAT  I  ON 

-  katz  e: 

< AMPLE 

RADIAN  FREQ. 

REAL  PART 

IMAGINARY  PART 

WGNITUDE 

PHASE(PAD) 

PHASE(OEG) 

.184761575 

-2.39228771 

-1.44991537 

2.91226675 

-2.53479993 

-145.227737 

.1149757 

-2.95134854 

-1.59799532 

2.54547622 

-2.59797991 

-143. 494148 

.124185488 

-1.76821733 

-1.34434233 

2.23461136 

-2.48373669 

-142.335942 

. 133488637 

-1.5329249 

-1.23744314 

1.97995693 

-2.46245977 

-141.333337 

.151991199 

-1.33732191 

-1.11928544 

1.74391144 

-2.44471727 

-148.372332 

.164819854 

-1.17441584 

-1.91985445 

1.54949324 

-2.43998665 

-139.235325 

. 183873829 

-1.93914292 

-.911228484 

1.38299273 

-2.42169891 

-133.753374 

.288923391 

-.924259954 

-.319533417 

1.23474443 

-2.41725149 

-123.493333 

.228513875 

-.331998732 

-.735817878 

1.11814154 

-2.41399932 

-133.541243 

.242812829 

-.7536395 

-.454941124 

.999379441 

-2.42424292 

-138.989983 

.24549878 

-.434823884 

-.584713337 

.992919543 

-2.43432243 

-129.591357 

.291585388 

-.438984245 

-.517732353 

.314122474 

-2.45443313 

-149.o2371 

.319924714 

-.533495553 

-.455338223 

.749989985 

-2.47382133 

-142.326851 

.351119174 

-.542791741 

-.397275798 

.472572198 

-2.59968973 

-143.794431 

.335352342 

-.397412444 

-.342944747 

.612443364 

-2.54721596 

-145.944774 

.422924291 

-.474353934 

-.292885857 

.553777222 

-2.59155843 

-143.434955 

.444158887 

-.44852322 

-.244322734 

.519755483 

-2.64231151 

-353.421999 

.589413884 

-.42394587 

-.199427917 

.44749443 

-2.79197383 

-154.769472 

.559881923 

-.399159742 

-.157212149 

.428995316 

-2.76638484 

-153.382233 

.413598733 

-.374138514 

-.117559888 

.394129531 

-2.33S78547 

-142.445991 

.473415872 

-.333595754 

-.9894242241 

.34242455 

-2.93795958 

-147.136313 

.739872211 

-.339992379 

-.9453343711 

.334941432 

-3.98395521 

-172.114937 

.81113984 

-.397735349 

-.8133945794 

.39S948847 

-3.89647223 

-377.414354 

.399215995 

-.233823344 

.9152296999 

.234234134 

3.83S91791 

174.939457 

.977889949 

-.259935971 

.9412779193 

.242393324 

2.93356371 

179.945954 

1.97224724 

-.233343953 

.3649954395 

.241942139 

2.37337962 

144.641232 

1.17481197 

-.29488473 

.3831222496 

.222953741 

2.75955234 

153.11374 

1.29154948 

-.179949323 

.3933763522 

.29597683 

2.44127174 

151.333777 

1.41747418 

-.152928919 

.199591363 

.183141559 

2.51979387 

144.373392 

1.55547417 

-.12437755 

.116795736 

.172921342 

2.39596163 

137.27S523 

1.79735247 

-.199884423 

.119814167 

.15643175 

2.27864179 

123.993238 

1.37331745 

-.8779413877 

.119135249 

.141923356 

2.14475823 

122.335439 

2.85451234 

-.3554543234 

. 1 15259292 

. 327936438 

2.81924576 

315.494391 

2.25791974 

-.934594795 

.133628593 

.114591277 

1.39593925 

198.577274 

2.4779744 

-.8294399183 

.3999475468 

.192924215 

1.77388134 

191.53553 

2.71353829 

-7. 588239 15E-93 

.3399511678 

.9992633153 

1.65393513 

94.7444937 

2.98364729 

2.34491874E-83 

.8793198393 

.8793514584 

1.53371942 

33.1423493 

3.27454922 

9.37769386E-93 

.3436214773 

.9493237529 

1.4273337 

31.383874 

3.59381373 

.9148346988 

.3533582695 

.9492147499 

1.32183415 

75.7354599 

3.94429ol3 

.9173137351 

.3438619345 

.3523895757 

1.22119313 

49.9443471 

4.32376137 

.8192341152 

.3493426436 

.9444931314 

1.12598582 

64.3394745 

4.75931926 

.8194712439 

.9328972762 

.938227753 

1.33434315 

59.3795424 

5.2148934 

.8183733774 

.926533782 

.9325414134 

.952529287 

54.575923 

5.72234773 

.8177244567 

.9211973958 

.3274324435 

.374333383 

59.3954394 

4.28829158 

.8162567436 

.9167939823 

.8233734391 

.891451719 

45.9332745 

4.39241226 

.8144344343 

.8132183633 

.9197159191 

.734297995 

42.9721911 

7.54443344 

.8129813711 

.3183277974 

.9145335321 

.472843491 

33.5952399 

3.39217537 

.8113794273 

3.33199572E-93 

.8139235434 

.41442121 

35.2157349 

9.11142777 

9.37985853E-83 

6.21376925E-93 

.8114741934 

.541735139 

32.1879273 

19.9888992 

3.51993595E-93 

4.79675333E-93 
Fig 

9.74963393E-83 
ure   20 

.5132213 

29.495425 
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Figure  20  (Cont . ) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

XX  BODE  PLOT  (AMPLITUDE)  XX 


ABSCI3SA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA 


=  . 1  RADIANS/SEC 
=  10  RADIANS/SEC 


AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-60  DECIEELS 


Figure  20 
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Figure  20  (Cont . ) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

XX  BODE  PLOT  (PHASE)  XX 


JL 


_L 


t   l_3— J I 


± 


JL 


_L 


I  I  I 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOM  MULTIPLES  OF  99  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  19  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  =  188  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  =  -189  DEGREES 


Figure  20 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  P ONERS  OF  S  OR  W 
-.  1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-.1  9 

1  9 

DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

9 

0 

4 . 6225 

4.3 

1 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 

-2.15  9 

-2.15  0 

0  0 

9  0 

MIN  GAIN  MAX  GAIN 

0  -30 

OPTION  TAKEN 

SIGMA  MIN  =  -1  SIGMA  MAX  =  1 

OMEGA  MIN  =  -1  CM EGA  MAX  =  1 


Figure  21   Compensated  System  Root  Locus 
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Figure  21  (Cont .  ) 


PROBLEM  IDENTIFICATION  -  KAT2  EXAMPLE 
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Figure  21 
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Figure  21  (Cont . ) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 
**  ROOT  LOCUS  PLOT  XX 
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ABSCISSA  ->  REAL  (SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHON  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -1  TO  2 

ORDINATE,  -1  TO  2 


Figure  21 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

This  thesis  demonstrates  that  common  classical  control 
programs  may  be  adapted  to  run  on  an  inexpensive  microcomputer 
system.   With  the  exception  of  the  memory  management  problem 
encountered,  this  was  done  with  relative  ease.   This  allows 
more  people  access  to  these  programs  by  introducing  a  new 
group  of  computers  on  which  they  can  be  run.   Also  these 
programs  require  no  knowledge  of  any  computer  language  or 
input  card  formats  making  them  easier  to  use  than  the 
Fortran  versions. 

It  is  also  demonstrated  that  the  transfer  function  pro- 
grams are  a  useful  tool  in  the  study  of  sampled  data  systems 
as  well  as  classical  systems. 

Only  a  sampling  of  five  programs  was  converted  in  this 
thesis.   Other  programs  can  and  should  be  similarly  converted 
to  run  on  microcomputer  systems.   There  are  programs  in 
existence  for  the  same  microcomputer  system  used  in  this 
thesis  that  generate  the  aircraft  stability  derivatives. 
An  effort  should  be  made  to  modify  these  programs  and  the 
programs  of  this  thesis  to  make  them  compatible  and 
complementary . 


164 


APPENDIX  A 
DESCRIPTION  OF  MICROCOMPUTER  SYSTEM 

The  microcomputer  system  used  in  developing  this  thesis 
consisted  of  the  following  components: 

Apple  II  plus  computer  (48K) 

Disk  II  5  1/4"  floppy  disk  drive  and  controller  card 

USI  9"  green  screen  monitor 

NEC  PC-8023A-C  dot  matrix  printer 

Grappler  printer  interface  card 

Add  Ram  16K  expansion  card 

The  programming  language  used  was  Applesoft  basic.   All 
graphs  were  generated  using  the  High  Resolution  graphics 
commands.   In  the  High  Resolution  graphics  mode  a 
matrix  of  dots  280  dots  wide  and  192  dots  high  can  be  dis- 
played.  The  High  Resolution  page  one  is  the  only  page  used 
in  this  thesis  and  it  resides  in  memory  in  the  8,192-byte  area 
from  $2000  to  $3FFF.   High  Resolution  page  two  resides  in 
the  area  from  $4000  to  $6000. 

A  memory  management  problem  was  encountered  during  the 
programming  of  the  programs  requiring  the  use  of  the  graphics 
capabilities  of  the  microcomputer  system.   To  understand  the 
nature  of  the  problem  encountered  a  brief  description  of  the 
normal  use  of  memory  in  the  Apple  II  is  necessary.   An  Apple- 
soft program  is  normally  loaded  at  memory  location  $800  and 
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loads  up.   LOMEM  is  set  to  the  end  of  the  program.   Simple 
variables  are  stored  from  LOMEM  upward  as  they  are  defined 
in  the  program.   This  gives  only  6K-bytes  of  RAM  before  there 
is  a  conflict  with  the  first  High-Resolution  page  and  only  14K 
before  there  is  a  conflict  with  the  second  High-Resolution 
page.   [Ref.  8]   Since  the  programs  and  variables  in  this 
thesis  exceed  14K,  simply  using  High-Resolution  page  two  in- 
stead of  one  is  not  the  answer.   The  first  part  of  the  solu- 
tion was  to  get  the  Applesoft  programs  to  load  above  the 
space  used  by  High-Resolution  page  one  and  in  effect  pro- 
tecting that  space  from  interference.   This  fix  creates  a 
new  problem.   The  disk  operating  system  loads  at  the  top  of 
the  48K  of  memory  and  sets  HIMEM  to  $9CF8.   String  variables 
start  at  HIMEM  and  build  down.   This  allows  22K  of  space  for 
the  program  and  all  variables.   Some  programs  in  this  thesis 
require  more  space  than  this.   Even  with  a  16K  memory 
expansion  card  installed  the  disk  operating  system  will 
ignore  it  so  the  added  memory  is  useless.   This  final  problem 
was  solved  by  using  a  utility  program  that  relocated  the  disk 
operating  system  into  the  higher  memory  provided  by  the 
memory  expansion  card  and  resets  HIMEM  to  $BFOO.   This  gives 
32. 5K  of  useable  program  space  which  is  sufficient  for  all 
programs  in  this  thesis.   [Ref.  9] 
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